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ABSTRACT 


The problem of finding the shortest path between two 
points in a plane containing obstacles is considered. The set 
of such paths is uncountably infinite, making an exhaustive 
search impossible. This difficulty is overcome by reducing 
the size of the search space. The search is first restricted 
to a countably infinite set by focusing attention on the set 
of homotopy classes. By applying simple optimality 
principles, a finite list of such classes is obtained whose 
union contains the shortest path. This process”) of 
simplification is accomplished by modeling the topology of the 
region with a graph. Optimality principles come into play 
during a graph traversal which is used to produce the finite 
search list. In addition, a computational investigation of 
two methods by which homotopy classes can be named is 
discussed, and properties of the graph models are 
investigated. The thesis of CPT André M. Cuerington, U.S. 
Army, calculates the actual shortest path using the search 


list produced here. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed 
in this research may not have been exercised for all cases of 
interest. While every effort has been made, within the time 
avallable, to ensure that the programs are free of 
computational and logic errors, they cannot be considered 
validated. Any application of these programs’ without 


additional verification is at the risk of the user. 
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i INTRODUCTION 


The industrial community relies today, in increasing 
measure, upon robots to handle a multitude of tasks. One goal 
is to develop robots, and fixed robot arms, with the ability 
to roam freely among obstacles. In this regard, the 


calculation of shortest paths is of obvious importance. 


A. THE PROBLEM 

As robot arms move, they must vary their configurations to 
position their end-effector or "hand". Once obstacles are 
introduced, the freedom of movement of the manipulator may 
become drastically reduced. The problem considered here is 
one of finding the shortest path from a starting point toa 
destination point in a planar region containing obstacles. 
This paper addresses several parts of the solution to this 


problem. 


B. OVERVIEW OF THE SOLUTION 

Between any two points in the plane, a path joining them 
may be chosen from an uncountably infinite set of 
alternatives. It is desirable to choose the shortest path 
from this set which, due to the presence of obstacles, may not 
be a straight line. An exhaustive search of the collection of 
possible paths is impossible, so another method must be 


developed. 


The set of possible paths is partitioned into a collection 
of equivalence classes--mutually exclusive subsets which 
collectively exhaust the partitioned set. This partition is 
induced when an equivalence relation is defined on the set of 
possible paths. This set of equivalence classes is countably 
Ine inate: 

Naming conventions are established which associate with 
each class a character string which allows the classes to be 
represented in a computer. Next, a finite list of candidate 
classes 1s produced to search. 

A heuristic is then applied to this finite list of classes 
to place them in an order which saves computational effort. 

The final step begins with this ordered list of candidate 
classes and solves the shortest path problem, class by class, 
in listed order. Savings of computational effort are realized 
when the class-by-class solution process can terminate without 
exhausting the ordered list of candidate classes. 

The key ideas in the above overview are now considered in 


somewhat greater detail. 


C. THE APPROACH 
l. The Topology 
Let P be the uncountably infinite set of all 
continuous, obstacle-avoiding paths from a starting point a to 
a destination point z. lee P; and p; are coterminal paths in 


P, that is, paths which have the same starting point and 


destination point, we say that Pp; 1s homotopic to p; if Pp; can 
be mapped to P. under a continuous function (with both 
endpoints fixed in place) without encroaching on any obstacles 
beet. lip. 223]. The homotopy relation is' reflexive, 
symmetric and transitive and therefore defines an equivalence 
relation [Ref. l:p. 223]. This relation induces a partition 
of P into a countably infinite collection of equivalence 
classes, known as homotopy classes [Ref. lip. 223]. 
2. Establishing the Naming Convention 

In order to name homotopy classes, a reference frame 
is established to represent the topological relationships in 
the plane with obstacles. For a given path, p, a string of 
characters, R(p), 1s recorded which encodes information 
concerning the relationship of the path to obstacles. 

Two algorithms are then presented which accept R(p) as 
input. These algorithms have the following important property 
which implicitly defines names for the homotopy classes: if 
p; and Pp, are coterminal paths and R(p;) and R(p;) are input to 
either of the two algorithms, then the outputs are identical 
if and only if P; 1s homotopic to P; [Ref. 2]. 

Appendix A presents a computational investigation 
which closely examines the methods employed by these two 
algorithms in producing the names for these homotopy classes. 

3. Generation of Names for Candidate Equivalence Classes 

An edge-labeled graph is constructed which models 


topological relationships within the region, where the nodes 


represent subregions induced by the reference frame. Nodes 
are connected if their corresponding subregions are adjacent. 
The names of the desired homotopy classes may be obtained by 
traversing this graph. This paper presents algorithms both to 
construct the graph and to traverse it. This traversal 
produces a list of candidate classes, which contain all paths 
of minimal length. 

The location of the origin affects the manner in which 
the plane is divided into wedges. In this regard, the graph 
which is created, given the location of the origin in one 


instance, may not be the same as the graph obtained if the 





Origin is moved to a new point. This thesis presents a method 
to count the number of possible different graphs, given n 
obstacles, and a method to construct them. 

4. Heuristic Ordering of Candidate Homotopy Classes 

For each class on the list, a lower bound on the 
length of its shortest representative path is constructed. 
The list of classes is then arranged into increasing order of 
these bounds. 

To obtain the bounds, a point is first chosen within 
each obstacle. A contraction deformation is then applied to 
each obstacle to “shrink"™ the obstacle so that the chosen 
point represents the obstacle. 

Then some class is fixed and its shortest path is 
examined. As all obstacles are simultaneously contracted to 


their representing points, this shortest path has a limiting 


position which is polygonal. The length of this polygonal 
path is the lower bound associated with that class. The 
length can be calculated from the class name and representing 
points without explrelt ly defining any contraction 
deformation. 
5. Class by Class Solution of the Shortest Path Problem 
In the final step of the solution, the classes on the 
ordered list of candidates are considered. The first class is 
removed from the list, l.e., that class with the smallest 
associated bound. The true shortest path in the named class 
is found by reversing the contraction previously applied to 
the obstacles, thereby transforming the polygonal path whose 
length provided the lower bound into the true shortest path. 
If the length of this path is smaller than the bound 
associated with the class on the top of the remaining list, 
the search is stopped. Otherwise, the first class is removed 
from that remaining list and the procedure is repeated. This 
procedure continues until the condition specified above is 


met, and the shortest path has been found. 


D. SUMMARY 

The solution to the problem of searching for the shortest 
path between two points in the plane with obstacles begins 
with consideration of a set of paths which is uncountably 
infinite. Through the homotopy relation, this set is reduced 


to a countably infinite set. The new set is further reduced 


to a finite list by modeling the regton with a Graphwand 
applying an optimality principle. This final list contains 
only those homotopy classes containing paths which are not 
self-crossing. From this list a shortest path is found. 
Computational effort is further reduced through the use of 
a heuristic which orders the list of candidate classes by 
increasing order of their lower bounds. The heuristic used 
also facilitates a methodical search of the classes while 
solving the shortest path problem. The use of this heuristic 
does not, however, imply that the solution 1s approximate. 


The solution to the problem will be exact using this method. 


E. THE CONTRIBUTION 

Chapter II of this thesis presents the methods used to 
establish a reference frame given the plane with obstacles and 
an algorithm which is used to generate the homotopy class name 
of a given path. 

Chapter III introduces a graph which is used to model the 
region. An algorithm is then presented which is used to 
traverse this graph to produce the list of candidate homotopy 
classes. Using the results of this algorithm, the thesis 
submitted by CPT André M. Cuerington, U.S. Army [Ref. 3], 
shows how to find the solution to the shortest path problem. 

Chapter IV introduces some properties of the graph which 


are used to model the region. 


Chapter V employs counting methods to determine the number 
of non-isomorphic graphs which may be encountered given a 
plane with some number of fixed obstacles. The number of non- 
1somorphic graphs represents the number of truly different 
reference frames which may be constructed in the region, thus 
producing different sets of names for the homotopy classes. 

Finally, Appendix A presents a computational investigation 
which addresses the question of whether or not two algorithms 
used to name homotopy classes actually produce the same 


results. 


Il. NAMING EQUIVALENCE CLASSES 


A. INTRODUCTION 

In this chapter, a labeling scheme is established to 
represent the equivalence classes. This notation is used 
throughout this paper to organize the computational search for 
the shortest path. 

An algorithm 1is~ presented that from a character 
representation of a path, p, names the path's respective 
equivalence class. It can be shown that, after being 
processed by the algorithm, different character 
representations of paths P; and P; yield the same output 
exactly when the two paths are in the same equivalence 
class (Ref. 2]. Thus, the name of a class will be the string 
obtained by applying the algorithm to any representative of 
the class. 

The computation of class names depends on a reference 
frame, Which in turn depends on a collection of obstacles. 
Although more than one reference frame can be drawn for a 
particular collection of obstacles, the choice of a particular 
frame fixes the representation of all homotopy classes. 

Once a reference frame is developed, the name of the 
homotopy class for a path can be determined by a two-step 
procedure. First, a character string R(p) is calculated which 


encodes certain information about the path taken through the 


obstacles. Second, this character string is accepted as input 
to an algorithm, which then produces a name for’ the 
equivalence class to which the particular path belongs. This 
algorithm produces the class names in terms of the same 
alphabet used to create the initial character string. The 
results from the algorithm described, Algorithm 1, will be 
used throughout this analysis as the class name associated 


with a given path. 


B. CONSTRUCTION OF REFERENCE FRAME 

Let b, be an arbitrary point in obstacle Bip = ee eer 
where n is the number of obstacles in the region. A point c 
1s chosen and a line drawn through each Dia tne ini Ee in extent 
in each direction, and having the following properties: there 
1S an open disk, &, centered at c such that 5NB, = ¢ for all 
k = 1,...,n, and the n lines connecting c with the points by 
are distinct. Such a point c can always be found as the above 
two conditions are satisfied by any point in the planar region 
that is neither on an obstacle nor on the n(n-1)/2 lines 
determined by pairwise choices of distinct b,. 


To draw a reference frame, n lines are first constructed 


joining c with each by. The line from c to by, is labeled as 
Ly . Each line is then partitioned into two directed, semi- 
infinite rays and one finite length line segment. The ray 


emanating from c in the direction away from by is called a,. 


The ray emanating from b, and away from c is denoted B, . The 


remaining line segment, [c,b, J, is also denoted a,. The 
reference frame 1s the collection of line segments and rays so 


constructed, as illustrated in Figure 2.1. 


B, b 
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Figure 2.1 A Reference Frame With n = 2 Obstacles 


C. RAW CHARACTER STRINGS 

A reference frame is fixed in a region T. Let a and z be 
points in T, and p be a directed path in T from a to z. Then 
the raw string of p, denoted R(p), is defined to be the 
ordered sequence of characters obtained by following p from a 
to z and recording, in order, the names of the rays that are 


crossed. 


Two special cases must be addressed to make the above 


definition complete. First, in the case that p crosses no 
reference rays, we set R(p) = e, where e denotes the empty 
Strange Second, if p crosses through c (simultaneously 
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the names of all a: will be recorded in order 


crossing all a, 


Siermcreasing, Subscript. 
Thus, raw strings are of the form 
Rigo) a Xp Xo ee Xp 


where each x; belongs to the alphabet 


A = ({e, Cece Ga, Bip. =, Bi). 
The character x; (Gaiam ee ee in eset neponamesco feathe jth 
reference ray crossed by p. Figure 2.2 shows a pair of paths 


connecting the two points a and z in the region with two 


obstacles. 





Figure 2.2 Paths a»B, and QQ 


D. ALGORITHM lL 
Algorithm 1 accepts as input a raw string R(p) for any 
path p in T and any reference frame. The output is a string, 


denoted C(R(p)), of characters also chosen from. the 


Tee 


alphabet A. This output string C(R(p)) @s thelname of Sta 
homotopy class to which p belongs. 

The algorithm is presented in terms of two functions. The 
first is the sorting function o. Let §S = X)-.-X, be any string 
over A. If S contains a two character substring Ke Xppp = iQy 
with i<j, then o(S) is the string which results by reversing 
the order of the leftmost such substring. Figure 2.3 depicts 
two such strings and makes clear that such paths are 


homotopic. If S contains no such two character substring 


then o(S) = S. 





esi 


Figure 233 Homotopie Al phawStemes QO} and QQ» 


Repeated application of o sorts all substrings’ that 
consist entirely of a,'s into non-decreasing order of 
subscript. E(S) is defined to be th@ String which results 


when this sorting is complete. 
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Hie  Secona EUmction is the cancellation function y. If 
string S contains a two character substring X, Xp4] with Kier 
then x(S) is the string which results by removing the leftmost 
such two character substring. Otherwise x(S) = S. Such 
cancellation 1s intuitive, for if a reference ray is crossed 
twice consecutively, then that is equivalent to not crossing 
at all. Thus repeated application of xy cancels all pairs of 
adjacent like occurences of identical characters. Let X(S) be 
defined as the string which results when all such possible 
cancellation is complete. 

With these definitions complete, Algorithm 1 is given in 


Figure 2.4. 


BEGIN 

k - 0 

5, o R(p) 

WHILE S, not equal to X(Z(S,)) 
k - ktl 
S) ~ X(£(S,.))) 

END WHILE 

Sk Gey =o; 

END 


Figure 2.4 Algorithm 1 
The output from Algorithm 1, C(R(p)), is called the 
Canonical representation of p and is the unique name for the 


homotopy class of p [Ref. 2]. That name is used throughout 


this paper. 
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I1I. A GRAPH MODEL FOR NAMING CLASSES 


This chapter presents a graph used to model the 
topological relationships in the region. The graph is then 
used to produce a finite list of names of candidate homotopy 
classes for subsequent searching. 


\' 


A. MODELING THE TOPOLOGICAL SPACE 

The reference frame, once established, divides the region 
into 2n "wedges". Each wedge is represented by a node. Two 
nodes are connected by one edge if the wedges they represent 
share a boundary of positive length corresponding to only an 
Q® ray, or two edges if the two wedges share a_ boundary 
corresponding to both a and 8B rays. No edge connects two 
nodes if the wedges ete represent have no common boundary 
other than the origin, c. Figure 3.1 illustrates this 
procedure of replacing the wedges by their associated graph 


Wise he no =e 





Figure 3.1 Establishing the Nodes and Edges of the Graph 
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Each edge of the graph is then labeled with an a, or a By . 
An a, 1s assigned to an edge if the incident nodes represent 
wedges sharing a boundary created by an ime eekerence line. 
The By labels are similarly assigned. When parallel edges 
connect adjacent nodes the outermost edge is labeled with By, 
the innermost edge with Oy - This labeling coincides with the 
fact that EJS Ibn Vey for each obstacle is on the outboard side 
of the obstacle with respect to the origin, c. Figure 3.2 


illustrates the edge labeling for n = 2. 





ool 
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Figure 3.2 Edge Labeling for Graph with n = 2 


In constructing the model, wedges and nodes are numbered 
counterclockwise using upper case Roman numerals. 

We define a ring graph to be such a graph used to model a 
region. Figure 3.3 illustrates the region with n = 2 and its 


associated ring graph. 


iS 





et 


Figure 3.3 Region and Ring Graph for n = 2 


If a and z are points in the topological space, T, ander 
is any path from a to z, then the string of characters, RG@poe 
Which represents that path can be produced by tracing p from 
a to z in T and recording R(p) as the list of rays crossed, as 
defined in Chapter II. However, R(p) can also be produced by 
traversing the ring graph and visiting the nodes’ which 
correspond to the wedges visited by p in T and recording the 
name of each edge as it is traversed. Figure 3.4 shows a path 
in a region where n = 2, and illustrates the path in the 
correspondine ring Gaarn. 

In constructing the model, it becomes unnecessary to 
generate a complete graph on the 2n nodes by eliminating the 
Orlgin as a boundary between any two wedges. This may be 
justified by observing that any path passing directly through 


the origin could be perturbed in such a manner that any 
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analysis of the path would remain unchanged. This situation 


momiblustrated ilmeFigure 3.5. 





Dil Oy 





R(p) = &@7 By &% 


Figure 3.4 Path Represented by Ring Graph for n = 2 





Figure 3.5 Perturbation of a Path Through the Origin 
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B. PRODUCING CLASS NAMES USING GRAPH TRAVERSAL 

Any path p in T has a character string representation 
associated with it which may be obtained by traversing the 
graph. A list of character strings which name the homotopy 
classes is desired. This list of names is finite and includes 
only the names of those classes whose representative of 
shortest length is not self-crossSing and starts and finishes 
in the appropriate wedges. 

A list of character strings is constructed by beginning at 
some node and searching the graph using a modified Breadth 
First Search (BFS) algorithm. This search utilizes a first- 
in/first-out queue which stores the current list of character 
strings and an element which denotes the node in which each 
string terminates. 

To illustrate the procedure, we begin with an empty string 
and the starting node of the graph on the queue. Then we copy 
the first string and the node denoting the current position 
from the queue. A new character which represents an edge 
which may be traversed in exiting the current node is added to 
the string. This new string is both printed and added to the 
bottom of the queue. All similar strings are generated by 
traversing each of the edges which may be traversed from the 
current node. The procedure is then repeated for the next 
string on the queue. Figure 3.6 presents this search in 


algorithmic language. 
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k := bottom of queue 

QO := top of queue 

v() := vertex on queue 
s() := string on queue 
BEGIN 

kta 

QO - ] 


v(k) + start vertex 
s(k) +e 
WHILE (Q<k) DO 
string - s(Q) 
vertex <- v(Q) 
FOR EACH EDGE OF THE FORM (vertex, vertex2) DO 
k «- k+l 
char + label[(vertex, vertex2)] 
s(k) «+ string:char 
v(k) + vertex2 
OUTPUT s(k), vk) 
Csr tL 
END WHILE 
END 


Figure 3.6 Modified BFS Algorithm 


In the absence of some criterion for terminating this 
search, an ever-growing list would be produced. A new string 
would be produced for every node which is adjacent to the 
current node. It is easily seen that this infinite list would 
eventually include the character string representation of 
every possible path in T which starts in some given wedge. In 
order to accomplish the task of reducing this list toa finite 
list of names of candidate homotopy classes, some stopping 
Criterion is required. 

By examining the list of strings, it can be seen that some 
OL the strings are not desirable because they are 


representative of classes which contain only paths which are 


nS 


self-crossing or are subject £6 cancel lation under 
Algorithm l. These paths should be removed from the list. 
Since all of the extensions of these excluded strings are also 
undesirable, they too should be excluded. The algorithm then 
terminates after producing a finite list. That the list is 
finite is substantiated by using the Pumping Lemma of finite 
state machines. This lemma states that for any language, an 
infinite string is possible if and only if cycles are 
permitted [Ref. 4:p. 77]. The language here consists of the 
R(p) strings. A cycle in the language is representative of a 
loop around any obstacle or set of obstacles. Therefore, if 
we remove the cycles, a language with a finite number of words 
1s produced, which in this case is a finite number of paths. 

Two such undeSirable strings are shown in Figure 3.7. 
These strings represent: paths which either wrap around an 
obstacle or cross themselves by circling two obstacles in 


figure eight fashion. 





Oy By a) a B) ay Boa) a» 


Figure 3.7 Paths Which Wrap and Cross 
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By performing several checks on each string as it is 
generated by the search algorithm, those strings which have an 
undesirable form may be "sieved" from the list. In other 
words, these strings will not be placed on the queue, and will 
therefore no longer be considered inthe search. The criteria 
which are applied to each string to accomplish this sieving 
are given below. 

The first criterion is that any string containing a 
substring of the form a, , By, ay, or 2, ,a,, 8, 1s not placed on the 
list. Such a string represents a class whose shortest 
representative wraps around an obstacle and is self crossing. 
The second criterion is that strings containing substrings of 
the form a, Bi, a), By or By, , a, B, , a, represent classes whose 
shortest representative is also self crossing. They too are 
not placed on the list. Figure 3.8 provides examples of paths 
which produce substrings of these types, and also demonstrate 
how the shortest representative in each class crosses itself. 

A third criterion eliminates from the list those strings 
which contain two like a's not separated by any B&B. Such 
strings are subject to cancellation as defined in Algorithm 1 
presented in Chapter II. Once cancellation takes place the 
resulting string merely duplicates some string already on the 
Prst. 

Another criterion introduces the restriction that the 
algorithm may not traverse an edge of the graph two times in 


immediate succession. Therefore, generation and cancellation 


ae 


of such strings 1s not required. Also, Gduplicatiensot stmine- 
CanenoEwOoceure 

Finally, those strings which do not terminate at the 
desired node are eliminated from the list. 

What remains on the list upon termination of the search 
given these sieving rules are the names of only those 
candidate homotopy classes which are to be used in the 


Subsequent. Seareh, for ushortest  patne 





B) 0B a) a) 


Figure 3.8 Examples of Self Crossing Representatives 
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IV. PROPERTIES OF THE RING GRAPH 


This chapter introduces some of the properties of the ring 
graphs obtained when the number of obstacles is greater than 
one. The case in which there 1s a single obstacle is 
addressed briefly at the end of the chapter. 

As stated in Chapter II, each reference line, Ly, in the 
region containing n fixed obstacles is partitioned into three 
parts--two semi-infinite rays and one finite line segment. 
Since there iS a one-to-one correspondence between these 
components and the edges of the ring graph, it 1s easily seen 
that the total number of edges in the graph will equal 3n for 
amy Th > 2. 

Now, we superimpose a ring graph over its corresponding 
reference frame, as illustrated in Figure 4.l. We first 
observe that, by definition, each Ly contains exactly one 
point b,. That point contributes to the existence of two 
edges connecting the adjacent wedges. Likewise, since no such 
I ities on L, on the opposite side of c, only one edge connects 
the adjacent wedges. Therefore, it is easily seen that in 
every ring graph with n 2 2, there are always two parallel 


edges opposite a single edge. 


ao 





Figure 4.1 Region and Ring Graph Superimposed for n = 3 


The number of edges incident with a node is called the 
degree of the node [Ref. 5:p. 86]. Consider the degree of the 
nodes of the ring graph. Every node ina given ring graph is 
of degree two, three or four. The sum of the degrees of 
the nodes ina graph is equal to twice the number of edges in 
that graph [Ref. 5:p. 86]. Since it has already been stated 
that the number of edges in the ring graph equals 3n, then for 


each ring graph, 


an 
degree(vertex i) = 6n. (4.1) 
1=1 
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Since we are interested in the degree of incidence of 


individual nodes we rewrite (4.1) in the form 


4 
6n = DE i(number of vertices of degree i) . (4.2) 
1=2 


It is an easy combinatorial result that in any graph, 
there are an even number of nodes of odd degree 
[Ref. 5:p. 91]. Therefore, in any ring graph, the number of 
nodes of degree three is even. 

Define g, h and j to be the number of nodes of degree two, 
three and four, respectively. Then the total number of nodes 
1s 


Geteie+ J = 2ne; (4.3) 


and (4.2) may be rewritten as 


29 +3h+47 = 6n. (4.4) 


Multiplying both sides of (4.3) by 3 and combining equations 


(4.3) and (4.4) yields 


3(9g +h+ j) = 29 + 3h + 47. (4.5) 


The solution to (4.5) shows that g = j or simply that the 
number of nodes of degree two equals the number of nodes of 


degree four in any ring graph with n 2 2. It may also be 
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verified by examining any ring graph that every node of degree 
two 1s directly opposite a node of degree four. Hence, there 
must be an equal number of each. 

All cycles in ring graphs are on an even number of nodes. 
Therefore these graphs may be classified as being of bipartite 
type. That is, the set of nodes may be partitioned into two 
sets such that all edges of the graph connect a node in one 
set to a node in the other set [Ref 5:p. 240]. Since the 
nodes are labeled sequentially in a counterclockwise 
direction, the two sets in this partition consist of those 
nodes whose labels are even integers and those whose labels 
are odd integers, respectively. 

This bipartite structure allows certain observations 
concerning the distance travelled in moving from node to node 
through the ring graph.’ Assigning an equal weight of one to 
each edge, the distance travelled along a path is simply the 
total number of edges traversed (with the multiplicity of 
traversals of individual edges included). It is easily seen 
that any path which begins and ends at nodes whose labels are 
of like parity will have even length. On the other hand, a 
path which begins and ends at nodes whose labels are of 
differing parity will have odd length. 

Consider a chain of m nodes of degree three. For m odd, 
the chain appears as that shown in Figure 4.2(a). For m even, 


the chain takes on one of the forms shown in Figure 4.2(b). 


Zo 


ee 
Bn ae 
Cp —a_ 0 


(a) (b) 

Figure 4.2 Chains of Nodes of Degree Three 
Consider a chain such as the one illustrated in 
Figure 4.2(a). It is easily seen that in concatenating this 
chain with nodes of degree two and four, it is possible to do 
so with only a degree two node on one end and only a degree 
four node on the other. Figure 4.3 illustrates the extension 


M@eesuch a chain. 


<a eo —c_ > — 


Figure 4.3 Concatenation of Even Degree Nodes 


Using a Similar construction, chains such as those 
illustrated in Figure 4.2(b) may be used to show, as seen in 
Figure 4.4, that only nodes of the same degree (two or four) 


may be concatenated on both ends. 
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Figure 4.4 Concatenation of Even Degree Nodes 


These facts lead to the observation that any string of 
nodes of degree two and any string of nodes of degree four 
must be separated by an odd number of nodes of degree three. 
Furthermore, if two nodes of degree two (or four) are 
separated by nodes of degree three, the number of degree three 
nodes must be even. 

It is possible to represent a given ring graph by an 
adjacency matrix, A. The resulting matrix 1s a symmetric 
2n x 2n matrix in which the (i,j) element is a lif there is 
a single edge connecting nodes v; and Vy, a 2 if two paral ven 
edges connect the nodes, or a0Oif the nodes are not adjacent. 
Figure 4.5 shows a ring graph with n = 3 and its associated 
adjacency matrix. By virtue of the conventions used to label 
the nodes, it is observed that all nonzero entries in the 
adjacency matrix are on the superdiagonal, the subdiagonal and 
in the (2n,1) and (1,2n) positions. (All entries on the main 
Giagonal of A are zero since no ring graph contains any edges 
which connect a node to itself.) Hence, the matrix becomes 


very sparse as n becomes large. 
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A basic property of an adjacency matrix is that the sum of 


th 


the values in the i” row (column) of the matrix is equal to 


Bie degree of node v. in the graph [Ret] Sipe se]. 

As stated above, the conventions used to label the nodes 
of the ring graph indicate that any two nodes whose labels are 
of like parity must be an even distance apart. By the 


bipartite property of the graph, all (1,3) entries in A, where 


1 and j are of like parity will necessarily be zero. 


a OO Qa 
2 C22 — a2 
ON Oe = 
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Figure 4.5 A Ring Graph and its Adjacency Matrix, A 


Another theorem from graph theory states that the number 
of paths of length m from node v; to node v, is the (1,j3) entry 
of the mt? power of the adjacency matrix, A" [Ref. 5:p. 144]. 
This theorem leads to the further conclusion that when m is 


odd, all (1,3) entries where i and j are of like parity will 


be zero. On the other hand, when mis even, all (i,j) entries 
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for i and j of opposite parity will be zero. Naturmatly, “agg 
entries not mentioned above may or may not be zero, and depend 
on the actual layout of the region and its corresponding ring 
graph. 

In the case where there is only one obstacle, n = l, a 
special case of the graph arises. Since one obstacle divides 
the region into two half planes, there are only two nodes in 
the graph with two parallel edges connecting them. Any two 
points in such a region may be connected by a path crossing 
either one a or one 8 edge. Therefore, this case is not 


considered in the analyses discussed in this paper. 
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V. NON-ISOMORPHIC RING GRAPHS 


A. INTRODUCTION 

The geometry of a region is fixed by the location of the 
obstacles. However, by varying the location of the origin, 
different reference frames may be established. These 
differences may lead to different ring graphs and therefore 
different homotopy class names. Figure 5.1 illustrates this 
possibility for n = 3, where two different reference frames 
result in different ring graphs and different class names for 
the same path. This leads to the following question: How 
many distinct ring graphs exist for a given collection of n 
fixed obstacles? 

It 1s shown in the sequel that the answer to the question 
may be found by determining the number of non-equivalent edge 
two-colorings of a 2n-gon subject to the restriction that 
Opposite edges of the polygon are colored differently. 
Burnside's Lemma is used to determine this number. 

It must be pointed out that coloring a graph in the 
context of this paper is not the traditionally accepted 
definition of a graph coloring. Here, adjacent edges (or 
later, vertices) may be colored the same color. 

Once the number of different ring graphs has been 
determined, Polya's Pattern Inventory Theorem is used to 


construct them. 


Si 


A complete example illustrating the procedures presented 
in this chapter is provided in Appendix C. Reference should 


be made to this example as new concepts are introduced. 
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Figure 5.1 Effect of Reference Frame on Ring 
Graph and Class Names 


B. THE RING GRAPH/POLYGON CONNECTION 
"A graph is completely determined by specifying its node 


and edge sets."[Ref. 6:p. 3] Any two graphs are isomorphic if 
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there exists a one-to-one correspondence between their node 
sets which preserves adjacency [Ref. 6:p. 4]. The question at 
hand then becomes: How many non-isomorphic ring graphs exist 
for a given collection of n fixed obstacles? 

Each ring graph associated with a problem containing n 
obstacles may be modeled using a 2n-gon whose vertices 
correspond to the nodes of the graph. The edges of the 
polygon are then colored as follows: a black edge (bold line) 
connects those vertices whose corresponding nodes are 
connected by two parallel edges; a white edge (fine line) 
connects those vertices whose corresponding nodes are 


connected by a single edge; and no edge connects those 


vertices which correspond to non-adjacent nodes. From this 
point forward, only single and double edges will be 
considered. 


It has been established that, in a ring graph, each set of 
parallel edges is opposite a single edge. Accordingly, the 
corresponding restriction on the polygon is that opposite 
edges must be colored differently. Figure 5.2 shows a ring 
graph for n = 2 and its associated 2n-gon. 

Since the configuration of a ring graph is defined by its 
node and edge sets, the location and/or spacing of its nodes 
is irrelevant. In constructing a graph then, the nodes may be 
thought of as being equally spaced around the ring. Therefore 
the 2n-gon associated with any ring graph may take the form of 


a regular polygon on 2n vertices. 
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Figure 5.24.Ring Graph andyZ2n-egon forenesa2 


With these ideas in mind, the problem of determining the 
number of non-isomorphic ring graphs given n fixed obstacles 
1s the same as determining the number of non-equivalent edge 
two-colorings of a regular 2n-gon subject to the restriction 


that opposite edges are differently colored. 


Caz COUNTING NON-EQUIVALENT EDGE TWO-COLORINGS 
l. Establishing the Permutations 

By coloring one edge of the 2n-gon, the restrictions 
which have been established previously automatically fix the 
color of the opposite edge. Therefore, of the 2n edges, n of 
them may be colored independently leading to a total of 2° 
possible colorings. 

Any pair of two-colorings of the polygon are 
considered equivalent if there is a rigid motion of the 2n-gon 


that maps the polygon onto itself {Ref 7:p. 3355 and wher 
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permutes one coloring into the other. Permissible rigid 
motions, or permutations, will be limited to rotations about 
the center of the polygon, reflections about diagonals drawn 
through opposite vertices, and reflections about bisectors 
drawn through opposite edges. Figure 5.3 illustrates each of 
these permutations using the square generated by n = 2. 
Examination of any 2n-gon reveals that there are a total of 2n 
rotations, n reflections about diagonals and n reflections 
about bisectors. If G is defined to be the set of all 
permissible permutations on the set of edges of the 2n-gon, 


then it follows that [G] = 4n. 
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Figure 5.3 Permissible Rigid Motions of a Polygon 


2. Fixed Colorings Under Given Permutations 
Any permutation may be factored into a product of 
G@isjoint cycles [Ref S:p. 691. In order for a coloring of the 
polygon to be preserved under that permutation, all edges in 


a cycle must have the same color. Since opposite edges are to 


85 


be colored differently, they must be in disjoint cycles. Then 
all edges which are opposite the edges of a given cycle must 
themselves be ina disjoint cycle, which will be called the 
companion cycle. In order for a coloring to be fixed under a 
given permutation, each cycle must have a unique companion 
cycle, hence the number of disjoint cycles must be even. In 
such cases the edges in one cycle are colored black while 
those in its companion'cycle are colored white. Figure 5.4 
provides an example of a coloring which is fixed under a 
permutation, and shows the companion cycles for that 


permutation. 
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Figure 5.4 Disjoint Cycles Produced by a Permutation 


We consider each of the three types of permutations. 
Throughout the following discussion ¥(x) is defined to be the 
number of edge two-colorings which remain fixed under a given 
permutation, n. For simplicity, the edges of each 2n-gon will 


be numbered clockwise beginning at some fixed edge. 
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a. Rotations 
The identity permutation, Ry, 1S that motion which 
does not move the polygon. Hence, ¥(n)) = 2". The remaining 
rotations are through the first 2n-1 multiples of (360/2n)° in 
the clockwise direction. 
Define a click to be (1/2n)'" of a full rotation 


of the polygon. RerOcactton, TM, 1 = 2...;C2n-L).. may be 


ie 
thought of as turning the polygon i clicks in the clockwise 
direction, thus moving each edge i positions along from its 
Starting pointe 

For the disjoint cycles generated by this 
permutation to have their coloring preserved, the number of 
disjoint cycles must be even. The edges in half of the cycles 
may be colored arbitrarily, while fixing the color of those in 
the companion cycles to the opposite color. Hence the number 
of choices to be made in coloring the edges of the polygon is 
one half the number of cycles, and ¥(n,) = gitcycles/2) | If the 
number of disjoint cycles is odd then (x; ) = 0. 

b. Reflections About Bisectors 

Once a bisector is drawn through two opposite 

edges of the polygon, the permutation for a reflection about 


that bisector yields a product of disjoint cycles consisting 


of two l-cycles and (n-1) 2-cycles. Such a permutation is 
illustrated in Figure 5.5. In order for each cycle to have a 
companion cycle, (n-1) must be even. Hence, n must be odd. 


In this case coloring half the 2-cycles fixes the color of the 


oy 


remaining companion cycles, and the number of possible choices 
for coloring the) (m-ljecycles ise@n— ey 2. One additional 


choice may be made in coloring those edges which lie on the 


bisector, resulting in a total of ¢c = (nt1)/2 choicesmio: 
coloring. Hence, when n is odd, ¥(n,) = 2°, where 
lim 2 gure nial 
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Figure 5.5 Disjoint Cyeles for Reflection About a Bisector 


For n even, the number of 2-cycles produced by the 


permutation is odd. This results ina cycle which fails to 
have a companion cycle. Hence no coloring may be preserved 
under the permutation. Thus, when n is even, P(x; ) = Q, 
le= 20 eon 1.) 


c. Reflections About Diagonals 
The permutation resulting froma reflection about 
an axis drawn through opposite vertices of the 2n-gon produces 
a product of disjoint cycles composed of nm 2-cyevece 


Figure 5.6 illustrates such ae reflection. It follows 
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immediately that when n is even, each cycle will have a 
companion and a total of n/2 choices may be made in coloring 
the 2n-gon. Then P(x; ) = 2nd i=! Spies An 

On the other hand, when nis odd, there will exist 
a cycle which does not have a companion cycle, and no coloring 
may be preserved under the permutation. This implies that 


when n is odd, » LD = Opa On ewer 1) 


. 
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Figure 5.6 Disjoint Cycles for Reflection About a Diagonal 


3. Burnside’s Lemma 
At this point, Burnside's Lemma is used to calculate 
the total number, N, of non-equivalent edge two-colorings of 
a 2n-gon given the restriction that opposite edges must be 
colored differently. The lemma is as follows. 
Let E be the set of edges of a regular polygon and let 
G be a group of permutations, XN, on the set E. Further, 
let C be any collection of colorings of E that is closed 


under G. Then the number, N, of non-equivalent colorings 
of the edges of the polygon is given by 
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e a 
ae Ta eck (ca 


where |G| is the number of permutations in G and ¥(x) is 
the number of colorings in C which are left fixed by 
Toker. 12 pemsas | 


Given the information derived earlier, this lemma may 


be invoked to produce N. This number is given by 
i 4n-1 
) Pix,) . (32) 
in b, POR 


Hence, the total number of non-isomorphic ring graphs given n 


fixed obstacles is obtained. 


D. PRODUCING NON-ISOMORPHIC RING GRAPHS 
Given the number of non-isomorphic ring graphs, WN, 
obtained by applying Burnside's Lemma, the next step is to 
Produce them. 
1. Using Partitions of the Integer n 
Givenn obstacles, we begin with a partition of n into 
summands, M,---,M,, where k is odd. Thus n = m,+m+...+m. 
With the summands in some fixed order, an edge two-coloring of 
the 2n-gon, and therefore a ring graph, may be produced in the 
following manner. Start at some fixed edge in the 2n-gon and 
GCOlor the fiimse m, edges black, proceeding clockwise from the 


starting edge. Color the next mM, edges white, the next m, 
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edges black, etc., until all k summands have been used. The 
colors of the remaining n edges become fixed as this procedure 
is followed. 

Using the procedure given above, only partitions with 
an odd number of summands are used. This guarantees that the 


th edge do not have the same color. 


starting edge and the 2n 
Therefore, the two-colorings of the 2n-gon are not equivalent, 
and the associated ring'graphs are not isomorphic. Figure 5.7 
illustrates two equivalent two-colorings of an octagon (under 
rotation through 45 degrees) using two different partitions. 


The first partition of n uses an odd number of summands, the 


second partition uses an even number of summands. 


Start Start 








Zee Ieee lterer lel 


Figure 5.7 Equivalent Two-colorings Resulting 
From Different Partitions 


It is possible that the number of candidate 


partitions, those with an odd number of summands, may actually 
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be less than N, the total number of non-isomorphic colorings 
calculated using Burnside's Lemma. When this occurs, there 
must exist at least one partition which, when ordered 
differently, is capable of producing more than one graph which 
is non-1isomorphic to the others. Figure 5.8 shows, forn=7, 
two non-equivalent two-colorings constructed from different 
orderings of the partition whose summands are 2,2,1,1,1. The 
question which arises" is: How many non-equivalent two- 
colorings, and hence non-isomorphic ring graphs may be 
produced from a given partition of n into an odd number, k, of 
summands? 

Polya's Pattern Inventory Theorem may be used to 


answer this question. 


Start Start 





i 2, Jae 


Figure 5.8 Non-Equivalent Colorings Constructed 
Using the Summands 2,2,1,1,1 of n = 7 
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2. Polya's Pattern Inventory Theorem 
An inventory is a list which indicates how many items 
®t a specific type exist [Ref. 7:p. 333]. This information 
may be encoded in the coefficients of a polynomial. Polya's 
Pattern Inventory Theorem follows. 
Let V be the set of vertices of a regular polygon. Let 
G be a group of permutations on the set V that acts to 
induce an equivalence relation on the colorings of V. The 


inventory of non-equivalent colorings of V using two 
colors 1s given by the generating function 


Pp ((b+w), (b’tw’),..., (b'twk)), 
where P, is the cycle index for the group of symmetries of 
the polygon. The inventory using colors Cie Cop ey 1s 
given by 
m m Mm 
PA (@mer me Cio, ec) 0s [Ref./ 7p. 355) — (5.3) 
j=l Jz=1 j=l 


3. Establishing the Cycle Index 

We begin with the set of all ordered partitions of n 
into k summands, where k is an odd integer. Any two 
partitions will be considered to be equivalent if there exists 
a cyclic shift of the elements of one partition, a complete 
reversal of the elements of one partition, or a combination of 
these which makes the partitions identical. Thus the set of 
ordered partitions of n may be separated into equivalence 
classes, with each class producing a non-equivalent two- 
coloring of the 2n-gon. 

The cyclic shifts and reversals of elements in a 


partition of k summands correspond to the rotations of a k-gon 
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and reflections of a k-gon about an axis drawn through a 
vertex and bisecting its opposite edge. Therefore, it is 
reasonable to consider the number of non-equivalent vertex 
colorings of a k-gon where the number of colors used is equal 
to the number of different summands in the partition. For 
example, with n = 7, a partition of n might be 2,2,1,7 0328 
Thus, k = 5. Here the vertices of a pentagon would be colored 
uSing two colors since only two different summands (1 and 2) 
appear in the partition. 

In examining any k-gon where k is odd, it is seen that 
k rotations and k reflections are possible for a total of 2k 
rigid motzons. The cycle index for the k-gon given these 
permutations is then produced. 

4. Applying Polya's Theorem 

Once the cycle: index, Pe, has been established, 
Polya's Theorem may be applied, where the number of colors, m, 
1S equal to the number of different summands in the given 
Partition. After the cycle index has been expanded using 
these m colors, the coefficient of the term whose exponents 
are equal to the multiplicities of each of the integers in the 
partition of n is determined. This coefficient represents the 
number of non-isomorphic ring graphs which may be produced 
from the given partition of n. 

The sum of these coefficients, after applying this 


procedure to all partitions of n with an odd number of 





summands, is equal to the total number, N, of non-isomorphic 
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ring graphs given n fixed obstacles. This represents the 
number of distinct reference frames which may be constructed 
in the plane with n obstacles, thus generating N different 


sets of homotopy class names. 


E. AN ALTERNATE SOLUTION TO THE COUNTING PROBLEM 

Upon completing the research and derivation of the 
solution method discussed above, an alternate solution was 
encountered in the literature. 

Fredricksen and Kessler [Ref. 9] present a procedure which 
may be employed in counting and generating the set of non- 
1somorphic ring graphs. This method is based upon an 
algorithm which produces the set of all lexicographic 


compositions of a positive integer, n. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Determination of the shortest path between two points in 
the plane containing obstacles has applications in many 
fields, particularly that of robotic path planning. This 
thesis presents a computational method which may be employed 
to generate the names of homotopy classes for a set of paths. 
By modeling the topological relationships of the region with 
a graph, a tool is made available for employing this method in 
producing only the names of those classes which contain 


candidates for the shortest path. 


B. RECOMMENDATIONS FOR FURTHER STUDY 

The field of graph theory is one that is gaining vast 
interest in mathematics and computer science. The graphs 
which are used to model the regions discussed in this thesis 
have several interesting characteristics. Further study and 
analysis of these graphs and their properties may aid future 
researchers in answering other questions concerning the 
shortest path problem. For instance, is there a sharp upper 
bound on the length of the class names which are generated by 
traversing the graph using the algorithm presented in 


Chapter <1 2 
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Another logical step in the shortest path problem is to 
extend the analyses presented here to include the third 
dimension. This extension follows since most industrial robot 
arms are of a pivoting arm type which are permitted to move 
freely through space while rotating at some fixed shoulder 
pont. 

The computational investigation presented in Appendix A 
was programmed with knowledge of only elementary FORTRAN 
programming skills. Although a complexity analysis was not 
done, it is believed that further study and refinement of the 
program code could result in a program which is much more 
efficient. 

Finally, this thesis introduces only part of the solution 
to the shortest path problem. The thesis submitted by CPT 
André M. Cuerington, U.S. Army, presents the remaining 
portions of the solution, and may be studied in conjunction 


with this paper. 
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APPENDIX A. A COMPUTATIONAL INVESTIGATION 


A. INTRODUCTION 

This appendix describes a computational investigation 
which was conducted as an aid in verifying that the homotopy 
class names produced by Algorithm 1, are reasonable when 
compared to a known method of naming homotopy classes 
[Ref. 2]. 

A second algorithm, Algorithm 2, is introduced here to 
forge a link between the class names of Chapter II and the 
well-known fundamental group. Two processes are presented to 
exploit this link: If Algorithm 1 truly names the classes, 
then both procedures should produce the same output for every 
path tested. One million test cases were examined and no 
counter-examples were found. While this computational 
evidence is not a proof, it does support the proof of the 


claim that Algorithm 1 names homotopy classes [Ref. 2]. 


B. ALGORITHM 2 
l. Fundamental Group 
The class name obtained from Algorithm 1 was written 
in terms of the alphabet A defined in Chapter II. We present 
an algorithm here where the name obtained is expressed in 


terms of the fundamental group of a topological space T. 
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The basic idea for producing a fundamental group for 
this problem is to regard the paths in Tas elements of the 
group, with path concatenation (*) as the group operation. 
There are two problems with considering the paths as the group 
elements. First, it 1s not necessarily possible to 
concatenate any two paths in T. In order to concatenate two 
paths using the operation (*), the first path must end at the 
point where the second begins. To make it possible to 
concatenate paths, in this analysis the test paths all begin 
and end at the point a. This point is called a base point. 

The second problem to overcome is that an inverse is 
not well defined for paths. If the identity is defined as the 
empty string, 1.e., stay at the base point, then when a path 
and its inverse are concatenated we have traced out a path 
which is not equal to the identity. To avoid this problem, 
the homotopy classes of paths are considered. Then all paths 
that have an empty character string representation will be in 
the identity class. Also, any class followed by its inverse 
will equal the identity. In Figure A.l, 9) represents the 
class a8, and (oi) RepLesenis Uthe hace Giaia 7, ce con) 
represents the class a) B; Bya, which reduces to the identity 
class by the rules of Algorithm 1. 

The group elements, therefore, are homotopy classes of 
loops around obstacles based at a common point, and the group 
operation is defined in terms of concatenation (*) of these 


classes. 
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Figure A.l1 A Homotopy Class and its Inverse 


Given a base point a, let p be any loop beginning and 
ending at a. Let [p;] denote the homotopy class of P;. So the 


set of group elements is ep. such |Ehae p; is a membenroma 


where P; is a loop based at a}. This set will be called G 
with the elements denoted g;- Figure A.2 illustrates several 
paths and their homotopy classes given n = 2 obstacles. 


Py 





Figure A.2 ([p)) = lp ers 


Now the group operation (*) can be defined by 


Ppl * ial = Sear 
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It is important to note that the fundamental group is 
finitely generated. As Figure A.3 illustrates, the generators 
are not unique. By fixing a set of generators, however, the 


class names become fixed. 





Figure A.3 Alternate Generators For Fundamental Group of the 
Space With Two Obstacles 


With the above information describing the fundamental 


group representation, Algorithm 2 is given below. 
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2. Functions in Algorithm 2 
a. Side Array 
Consider the reference lines L, constructed jin 
Chapter II to be oriented rays in the direction from c to by . 
It then becomes reasonable to discuss a ‘right' anda ‘left' 
side of those lines. We allow a moving point to trace a path 
p from start point to destination point. Side(k) is a 
function which defines the side of L, on which the moving 
point lies. The output is either ‘left' or ‘right'. The 
output is never ‘on’ because this routine is used only after 
a complete crossing of Ly takes place. A crossing is 
considered to be complete when the moving point leaves L, to 
one side after having met Ly from the opposite side. 
b. Switch Function 
The switch function is defined by: 
( “right” if side (k) =" tet te 
switch (side(k)) = 4 
| ‘Yeft" if side (kK) = %richius 
As the moving point traces p, each time a reference linesig 1s 
crossed, the switch function will be applied to indicate on 
which side of L, the moving point lies. 
c. Index Function 
Let Xy be any character from the alphabet A 
representing the jth element in R(p). Let index(x) = j if x 
equals a: or B.; thus, the jth crossing completed by the moving 


j 


point 18S a crossing of L So if the moving point is on the 


j e 
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left side of L, initially and this line is now crossed on 
either side vot b,7 then side(3) would equal ‘right’. 
d. The Algorithm 

Using the functions described above, Algorithm 2 
is given in Figure A.4. The algorithm works in two phases. 
The first phase initializes the array side(k) for k=l,...,n. 
The second phase reads the raw string from left to right and 
adds an element to the fundamental group representation for 
each 8B crossing. When the end of the input string is reached 
the output, F(R(p)), is a shorter string with one character 
corresponding to every 8 element in the original raw string. 


Each character represents a generator or its inverse. 


begin 
inpure (Dp) = XX. +. + Xp 
MOY j=! yee 


If ais to right of reference line L: then 
J 
side (j) «+ right 
else 
side (j) - left 
end if 
end for 
CG - ja] 
for Kee 1.4. mM 
1a cex cy) 
side(i) + switch(side(i) ) 
1f x, = 8, for some r, then 
if side(index(x,)) = left, then 
Ge-G* J; 
else 
G+ 6G * (g,)7} 
end if 
end if 
end for 
WGEGE) ees) G 
end 


Figure A.4 Algorithm 2 


ao 


e. Fundamental Group Cancellation Function 

The cancellation rules in the fundamental group 
differ somewhat from those for the raw string. Although sets 
of generators are not unique, we can obtain a_e unique 
representation of each class with respect to a particular set 
of generators. For any given set of generators, every 
homotopy class can be represented as a product of these 
generators and their inverses. Even so, this representation 
is not unique until cancellation 1s applied. The cancellation 
rule follows. 

Let G= {Q),---, gn} be a set of generators of the 
fundamental group of T (base point a) and let Y = Yi Yo. - wy be 
a representation of some homotopy class in terms of the g; and 
their inverses ™ 1-em@ior jeach 1 - 1) ny. Q; or (chee 
Lor some 3, kK 2-52). 

The cancellation function Kx is defined as follows. 
If Y contains a two character substring y:yi,) with iG (Ye 
then K(Y) is the string which results by removing both y; and 
Yay in their leftmost occurrence, from Y. Finitely many 
repeated applications of kK produce a string in which no 
further cancellation 1s possible. We define this string as 


K(Y). 
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C. THE COMPUTATIONAL INVESTIGATION 
1. The Approach 

A graphical representation of the two procedures is 
provided in Figure A.5. The investigation begins by accepting 
a random path p as input and generating two names for the 
homotopy class representing p which are output. 

The idea of the test is to determine for a given path 
its raw string character representation. For computational 
purposes, the paths considered are polygonal paths. This raw 
string is then input to both algorithms. 

The raw string (R(p)) is first entered into 
Algorithm 2. The output (F(R(p))) from Algorithm 2 is then 
subject to the cancellation routine which outputs the fully 
cancelled string (K(F(R(p)))). Next, R(p) is input to 
Algorithm 1 which produces the canonical string (C(R(p))). 
fepgorithm 2 is then applied to C{R(p)) to output F(C(R(p))). 
The results of these tests are then compared, and the pogman 
continues to run until a predetermined number of paths are 
checked. 

It can be shown that C(R(p)) is a unique class name if 
and only if F(C(R(p))) equals K(F(R(p))) for all p in 
T [Ref. 2]. Before the proof of this conjecture was obtained, 
the above procedures were programmed and tested for one 
million different paths and no counter-examples were 


discovered. 


D) 


R(p) 
ALG | AiG e7 
Cena is F(R(p)) 
AG 7 K(FCR(p))) 
meCCRi p))) K(F(R(p))) 


Figure A.5 Flow Chart of Two Algorithms 


2. The Test 

The general algorithm used to test the two procedures 
represented in Figure A.5 is easily followed ane is shown in 
Figure A.6. However, the attached program, which was used to 
test the model, is substantially more involved. This was done 
in an attempt to write efficient programming code. 

Two time saving techniques are employed in order to 
test the one million cases in 25 minutes. The first method 
involves subroutines which perform cancellation. It is simple 
to code a program that scans the character representation 
repeatedly to find all possible cancellation in the string 


BB, a) a, By, which reduces to B,. However, in the enclosed code, 


ofS 


pointers are used to mark the position where the first 
cancellation occurs and then check the newly adjacent 
characters for cancellation. This method reduces’ the 
complexity of an algorithm and is applied in both Algorithm l 


and Algorithm 2. 


mn = number of reference frames to be considered 
m = number of paths to be considered 
nobs = number of obstacles on each board 
seed = input seed for random generator (not needed?) 
begin 
fOr boakra., = 1o4-.4. 1 
create reference frame or board 
PE OrmoOdat ha =arl , . ..04.1 2 
1. Create a polygonal path p with nseg 
segments 
2 Form R(p) 
3. Form F(R(p)) 
4. Form K(F(R(p))) 
5 Form C(R(p) ) 
6. Form F(C(R(p))) 
if K(F(R(p))) does not equal F(C(R(p))) then 
print (raw string, board, and path info 
for the counter example) 
end if 
end for 
end for 
print (final seed) 
end 


Figure A.6 Test Algorithm 


The second time-saving device involves the sorting of 


2 


character strings. Instead of a bubble sort method (order n 
complexity), a merge sort algorithm (order n(log(n)) 
complexity) is used. As can be seen in the attached code, 


this choice to save computer time calls for a significant 


i) 


increase in programmer time, which is generally much more 
expensive. 
The code which is used to implement this test is 


presented in the remaining pages of this appendix. 
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000 anN9ad000 


9N09N0 0000 


PROGRAM CLASNAME FORTRAN 


ROR AI ORR OR A OR COR RR OO IRR RR AOR BR OR AO OR FOR A TOR OOK ARR ARR A 
THIS PROGRAM RANDOMLY GENERATES OBSTACLES AND A POLYGONAL PATH 
THROUGH THOSE OBSTACLES TO DETERMINE WHETHER OR NOT THE PROCEDURES 
USED BY COMPETING ALGORITHMS PRODUCE THE SAME FUNDAMENTAL GROUP 


REPRESENTATION FOR THAT PATH. 
HAE EMILIO BE A IR YOM CRE CAE SK aK PEE A IC RC LO IR A AEB AR A IC oR ICR A IE CIR RE A CACO I OE FOI Co 


THIS PORTION OF THE PROGRAM SERVES AS THE MAIN DRIVER WHICH RECEIVES 
THE PARAMETERS DEFINING THE REGION AND INITIALIZES THE LINK LIST 
ARRAYS WHICH WILL BE USED TO REPRESENT THE POLYGONAL PATH. 


INPUT: INITIAL SEED FOR THE RANDOM NUMBER GENERATOR, NUMBER OF 
OBSTACLE CONFIGURATIONS AND PATHS TO BE TESTED, NUMBER OF 
OBSTACLES IN THE PLANE AND NUMBER OF SEGMENTS IN EACH 
POLYGONAL PATH 


OUTPUT: FINAL RANDOM NUMBER GENERATOR SEED AND MESSAGES INDICATING 
ANY PATHS WHICH PRODUCE DIFFERENT FUNDAMENTAL GROUP 
REPRESENTAT IONS 


REAL*8 BX(1000), BY(1000), X(1000), Y( 1000), DSEED 
INTEGER NOBS, NSEGS, NUMPTS, N, M, HEAD( 1000) 

+ NEXT( 1000), PRED( 1000), BOARD, PATH 
EXTERNAL GGUBFS 


DSEED = 123457.0 


NSEGS = 5 
= NSEGS + 1 
PRINT*,‘ INPUT SEED‘, DSEED 


DO 1 BOARD = 1, N 


CALL BOARDS (DSEED, BX, BY, NOBS) 
DO 2 PATH = 1, M 


CALL INIT (HEAD, NEXT, PRED, PSEED, DSEED ) 
CALL PATHS (DSEED, X, Y, NSEGS) 
CALL TEST(HEAD,NEXT,PRED,BX,BY,NOBS,X,Y,NUMPTS , BSEED,PSEED) 


CONT I NUE 
OONT | NUE 


— BO) 


PRINT¥," ' 
PRINT*, ‘FINAL SEED', DSEED 


STOP 
END 
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OO 909 20 


SUBROUTINE INIT ( HEAD, NEXT, PRED , PSEED, DSEED) 
RA RMI I FR A RIOR 2 MR AR IT IR A A RE I IR CR AOR IE OR A RE ORE AOR RAK RRR 


THIS SUBROUTINE INITIALIZES THE DOUBLE LINK LIST ARRAYS THAT WILL 
REPRESENT THE PATH. 


INPUT: DSEED 


OUTPUT: HEAD, NEXT, AND PRED ARRAYS SET TO ZERO, AND SEED FOR THE 
R.N.G@. PRIOR TO CONSTRUCTION OF THE PATH 


INTEGER HEAD( 1000), NEXT( 1000), PRED( 1000) 
REAL*8 PSEED, DSEED 


PSEED = DSEED 


DO 3 K = 1, 1000 
HEAD (K ) 
NEXT (K) 
PRED(K ) 

3 CONTINUE 


ee 


0 
0 
0 


RETURN 
END 


SUBROUTINE BOARDS (DSEED, BX, BY, NOBS) 
HAA H HA HH HA HK HH HH HH AH HH HH HK HK HH HEH HK A HK KH HK HK IH HHH HAC HK HE HH KC I HK HK IK A Hk HK IK 


THIS SUBROUTINE USES A PSEUDO RANDOM NUMBER GENERATOR TO CREATE THE 
COORDINATES OF EACH OBSTACLE ON THE BOARD, SCALING ALL COORDINATES 
TO BE IN THE INTERVAL (-1,1). 
INPUT: NUMBER OF OBSTACLES IN THE REGION AND A SEED FOR THE R.N.G. 
OUTPUT: OBSTACLE COORDINATES 


REAL*8 BxX(1000), BY( 1000), DSEED 
EXTERNAL GGUBFS 


DO 1 1 = 1, NOBS 
BX(1) = 2. * GGUBFS(DSEED) - 1. 
1 BY(1) = 2. * GQUBFS(DSEED) - 1. 
RETURN 
END 


SUBROUTINE PATHS (DSEED, X, Y, NSEGS) 


CFR RRA ARE A a RR a A RR OG Ro aR AO RR A a AR ROK a ORK RK RR ORR ORO I aa RR OR OR RR 


Oo O98 9090900 


Ord "OOO © Bo 


© 


RANDOMLY GENERATES THE X AND Y COORDINATES OF THE VERTICES FOR THE 
POLYGONAL PATH, SCALING ALL COORDINATES TO BE IN THE INTERVAL (-1,1). 
THIS SUBROUTINE ALSO ENSURES THAT THE PATH |S A CLOSED LOOP BY 
ASSIGNING THE START/FINISH POINTS THE SAME COORDINATES. 


INPUT: NUMBER OF PATH SEGMENTS AND A SEED FOR THE R.N.G. 
OUTPUT: COORDINATES OF VERTICES ALONG THE POLYGONAL PATH 


REAL*8 X(1000), Y( 1000), DSEED 
EXTERNAL GGUBFS 


DO 11 = 1, NSEGS 
X(1) = 2. * GQUBFS(DSEED) - 1. 
1 Y(1) = 2. * GGUBFS(DSEED) - 1. 
X(NSEGS+1) = X(1) 
Y(NSEGS+1) = Y(1) 
RETURN 
END 


SUBROUTINE TEST(HEAD, NEXT, PRED, BX, BY, NOBS, X, Y, NUMPTS, 
+ BSEED, PSEED) 


TRAE RARE PE EEE SESE SES SEE OE PCT oR COSI IR OR SR OR OKC aK HR aC CIR 2 CC IC a OR OR IC OC RCE FE I 2, oR ap oe 


GENERATES THE RAW STRING OF CHARACTERS REPRESENTING A PATH AND 
DETERMINES ITS FUNDAMENTAL GROUP REPRESENTATION USING ALGORITHMS 1 
AND 2. RESULTS OF THESE COMPETING ALGORITHMS ARE THEN COMPARED FOR 
DIFFERENCES. 


INPUT: A POLYGONAL PATH AND COORDINATES OF POINTS REPRESENTING 
OBSTACLES 


OUTPUT: COMPARED RESULTS OF ALGORITHMS 


INTEGER HEAD( 1000), NEXT(1000), PRED( 1000), FR( 1000), NR( 1000), 
+ PR( 1000), FL( 1000), NL( 1000), PL(1000) 
REAL*8 8xX(1000), BY(1000), X(1000), Y( 1000) 


DATA FR/1000*0/ 
DATA NR/1000*0/ 
DATA PR/1000*0/ 
DATA FL/1000*0/ 
DATA NL/1000*0/ 
DATA PL/1000*0/ 
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CALL RAWSTR(BX, BY, NOBS, X, Y, NUMPTS, HEAD, NEXT, PRED, NELEMS) 
CALL ALG2 (NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FR, NR, PR) 
CALL CANALG2 (FR, NR, PR) 

CALL ALG1 (HEAD, NEXT, PRED, NELEMS) 

CALL ALG2 (NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FL, NL, PL) 
CALL CTREX (FL, NL, PL, FR, NR, PR, BSEED, PSEED) 


RETURN 
END 


SUBROUTINE RAWSTR(BX, BY, NOBS, X, Y, NUMPTS, HEAD, NEXT, PRED, 
+ NELEMS ) 
EA KK KK KK EE KEK EK KK EK KK KK IK EE KK KK KK KK KK HK KK RK HEM HK HK KKK KK KEI IK KK 3K 9K HK KK IKK IK 


PRODUCES THE RAW STRING OF CHARACTERS WHICH REPRESENTS THE PATH, 
PAYING CLOSE ATTENTION TO THE ORDERING OF THE CHARACTERS WHERE 
NECESSARY. THE STRING IS CONSTRUCTED BY IDENTIFYING THOSE OBSTACLE 
REFERENCE LINES WHICH ARE CROSSED AS EACH PATH SEGMENT IS TRAVERSED 
IN ORDER. 


INPUT: ORDERED LIST OF VERTICES REPRESENTING THE PATH AND A POINT 
REPRESENTING EACH OBSTACLE IN THE REGION 


OUTPUT: A RAW STRING OF ALPHAS AND BETAS CONTAINED IN AN ARRAY 
NAMED HEAD() AND ITS PARALLEL ARRAYS NEXT() AND PRED() WHICH 
PRODUCE THE DOUBLE LINKED LIST 


IMPLICIT REAL¥8 (A-H,0O-Z) 

INTEGER HEAD, NEXT, SEG, FSTSTR, HOINDX, SEGEND, FIRST, START, 
+ LENGTH, PRED 

LOGICAL ALLALF 

DIMENSION X( 1000), Y( 1000), Bx(1000), BY( 1000), A1(1000), 
+ B1( 1000), A2(1000), B2( 1000), D2( 1000), DIST( 1000), 
+ FIRST( 1000), HEAD( 1000), NEXT( 1000), PRED( 1000) 


DATA DIST/1000*0 .0/ 
DATA FIRST/1000*0/ 


CALL SETUP (NUMPTS ,NOBS ,CX,CY,X,Y,BX,BY,A1,B1,A2,B2,D2) 


HDINOX = 2 
NSEGS = NUMPTS - 1 
LOLD = 1 


HEAD(LOLD) = 0 
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DO 7 SEG = 1, NSEGS 


SEGEND = SEG + 1 
FSTSTR = HOINOX 

LSTSTR = HOINOX 

ALLALF = .TRUE. 

LENGTH = 0 


DO 6 LINE = 1, NOBS 
CHECK1 = (A1(LINE)*X(SEG)+B1(LINE)*#Y(SEG) )* 
+ (A1(LINE )*X(SEGEND) +B 1(L INE )*Y(SEGEND) ) 
IF (CHECK1.LT.0) THEN 
LENGTH = LENGTH + 1 
CHECK2 = (A2(SEG)*BX(LINE)+B2(SEG)* 
+ BY(LINE )+02(SEG) )*D02( SEG) 
IF (CHECK2.LT.0) THEN 
HEAD(HDINDX) = -LINE 
LSTSTR = HDINDX 
HDINOX = HDINDX + 1 
ELSE 
CALL CASES1 (A1, B1, A2, B2, 02, SEG, LINE, XINT, 
+ YINT) 
DISTC = XINT¥¥2 + YINT¥*2 
DISTB = (XINT-BX(LINE) )¥**2 + (YINT-BY(LINE) )¥**2 
IF (DISTB.LT.OISTC) THEN 
HEAD(HDINDX) = LINE 


LSTSTR = HDINDX 

HDINDX = HDINDX + 1 

ALLALF = .FALSE. 
ELSE 


HEAD(HDINOX) = -LINE 
LSTSTR = HDINDX 
HDINDX = HDINDX + 1 
ENDIF 
ENDIF 
ENDIF 
6 CONT | NUE 
IF (LENGTH.NE.O) THEN 
IF (ALLALF) THEN 
CALL ALPHAS (NEXT, LOLD, FSTSTR, LSTSTR, SEG, START, 
+ HEAD, NSEGS, PRED) 
ELSE 
CALL ORDER(SEG, LINE, HDINDX, HEAD, NEXT, Al, A2, 
r B1, B2, 02, FSTSTR, LSTSTR, X, Y, START, 
r LENGTH, FIRST, LOLD, NSEGS, PRED) 
ENDIF 
ENDIF 
7 CONTINUE 


CALL COUNTR(START, NEXT, HEAD, NELEMS) 


RETURN 
END 
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SUBROUTINE SETUP(NUMPTS, NOBS, CX, CY, X, Y, BX, BY, Al, Bl, A2, 
+ B2, D2) 


POR HHT RR EO A FE RE RRR ROKR ER RRR RRA ROKR 


FOR EACH OBSTACLE, THIS ROUTINE DETERMINES THE COEFFICIENTS OF THE 
FQUATION FOR THE REFERENCE LINE FROM THE OBSTACLE TO THE ORIGIN. IN 
ADDITION, IT CALCULATES THE COEFFICIENTS OF THE LINE REPRESENTING 
EACH SEGMENT OF THE POLYGONAL PATH. 


INPUT: NUMPTS, NOBS, COORDINATES OF VERTICES ALONG POLYGONAL PATH 
AND COORDINATES OF THOSE POINTS WHICH REPRESENT EACH 
OBSTACLE 


OUTPUT: COEFFICIENTS OF LINEAR EQUATIONS REPRESENTING PATH SEGMENTS 
AND REFERENCE LINES FOR EACH OBSTACLE 


IMPLICIT REAL*¥8 (A-H, O-Z) 
DIMENSION X(1000), Y( 1000), BX(1000), BY( 1000), A1i(1000), 


+ B1( 1000), A2( 1000), B2(1000), D2( 1000) 
DO 4 | = 1, NOBS 
Al( |) = BY(I) 
B1CI) = -BX(1) 
4 CONTINUE 
CX = 0.0 
CY = O20 


4 


DO 5 J = 2,NUMPTS 
A2(1) =ay(J) — YC) 


B2(1) = Xl) =m) 
D2(1) = CYCI)*XCI) )-(XC1)FYCJ)) 
| = 1+1 
5 CONTINUE 
RETURN 


END 
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SUBROUTINE CASES1(A1, B1, A2, B2, D2, SEG, LINE, XINT, YINT) 
KMART RRA EK I HE FR IA RR OR IR EAR A OR RR OR ORR 


THIS SUBROUTINE DETERMINES THE COORDINATES FOR THE POINT OF 
INTERSECTION OF A GIVEN PATH SEGMENT AND A GIVEN OBSTACLE REFERENCE 
LINE. NUMERICAL STABILITY OF CALCULATIONS REQUIRES THE MATHEMAT! CAL 
OPERATIONS TO BE SEPARATED INTO CASES, THE BEST CASE BEING USED FOR 
EACH PARTICULAR SITUATION. 


INPUT: COEFFICIENTS OF LINEAR EQUATIONS FOR PATH SEGMENT AND 
REFERENCE LINE TO BE EVALUATED 


OUTPUT: POINT OF INTERSECTION OF THE PATH SEGMENT AND THE OBSTACLE 
REFERENCE LINE 


IMPLICIT REAL¥8 (A-H,O-Z) 
INTEGER SEG 
DIMENSION A1(1000), B1(1000), A2( 1000), B2(1000), 02( 1000), A(2,2) 


A(1,1) = A1(LINE) 
A(1,2) = B1(LINE) 
A(2,1) = A2(SEG) 
A(2,2) = B2(SEG) 
BIGEST = 0.0 
Domdet = 1,2 
Dom 6 la 
TEST = DABS(A(K,L)) 
IF (TEST.GT.BIGEST) THEN 
BIGEST = TEST 
KBIG = K 
LBIG =L 
ENDIF 
1 CONTINUE 


IF (KBIG.EQ.1) THEN 
IF (LBIG.EQ.1) THEN 
YINT = -D2(SEG)/(B2(SEG)-B1(LINE)*A2(SEG)/A1(LINE) ) 
XINT = -B1(LINE)¥*YINT/A1(LINE) 
RETURN 
ELSE 
XINT 
YINT 
RETURN 
ENOIF 
ELSE 
IF (LBIG.EQ.1) THEN 
YINT = (02(SEG)*A1(LINE)/A2(SEG) )/ 
+ (B1(LINE)-B2(SEG)*A1(LINE)/A2(SEG) ) 
XINT = (-D2(SEG)-B2(SEG)*Y INT) /A2(SEG) 
RETURN 


-D2(SEG)/(A2(SEG)-A1(LINE)*B2(SEG)/B1(LINE)) 
-A1(LINE )*X1INT/B1(LINE) 
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ELSE 
XINT = (D2(SEG)*B1( LINE) /B2(SEG) )/ 


‘ (A1(LINE )-A2(SEG)*B1(LINE )/B2(SEG) ) 
YINT = (-D2(SEG)-XINT*A2(SEG) )/B2(SEG) 
RETURN 
ENDIF 
ENDIF 
END 


SUBROUTINE ALPHAS (NEXT, LOLD, FSTSTR, LSTSTR, SEG, START, HEAD, 
+ NSEGS, PRED) 
apo coke) fark anos eS So oo eink ok gk ok Ge ao kok ok ac ak ek ok ie aka ak 


GIVEN THAT A SEGMENT OF THE PATH CROSSES ONLY ALPHA RAYS WHEN 
TRAVERSED, THESE CROSSINGS ARE SIMPLY INSERTED INTO THE LINK LIST 
ARRAYS IN THE ORDER IN WHICH THEY WERE DETECTED, |.E., SMALLEST TO 
LARGEST IN ABSOLUTE VALUE. 


INPUT: STRING REPRESENTING THE REFERENCE LINES CROSSED BY THE 
CURRENT PATH SEGMENT, GIVEN THAT ALL CROSSINGS ARE ALPHAS 


OUTPUT: UPDATED NEXT() AND PRED() ARRAYS WHICH CONTAIN THE STRING 
REPRESENTING THE MOST CURRENT PATH SEGMENT 


INTEGER NEXT( 1000), FSTSTR, SEG, START, HEAD(1000), PRED( 1000) 
NEXT(LOLD) = FSTSTR 

PRED(FSTSTR) = LOLD 

LAST = LSTSTR - 1 


IF (LAST.GT.FSTSTR) THEN 


DO 12 | = FSTSTR,LAST 
NEXT(!1) = | + 1 
PRED(I+1) = | 
12 OONT | NUE 


ELSEIF (LAST.EQ.FSTSTR) THEN 
NEXT(FSTSTR) = LSTSTR 
PRED(LSTSTR) = FSTSTR 

ENDIF 


IF(SEG.EQ.1) START = LOLD 
LOLD = LSTSTR 


RETURN 
END 
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SUBROUTINE ORDER(SEG, LINE, HDINDX, HEAD, NEXT, Al, A2, Bl, B2, 
+ DZ, FSTSTR, LSTSTR, X, Y, START, LENGTH, FIRST, 
a LOLD, NSEGS, PRED) 


SRR ICIORIOIIOI ECR IOI ain ICICI ICR CR aR ak a RIOR IC aK ac ok a6 a ak I ka a eS icc 3c acai Si age a aa 


GIVEN THAT A SEGMENT OF THE PATH CROSSES ONE OR MORE BETAS, THIS 
SUBROUTINE DETERMINES THE ACTUAL ORDER OF CROSSING WHEN TRAVERS!NG 
THE SEGMENT FROM BEGINNING TO END. THIS IS DONE BY FIRST DETERMINING 
THE DISTANCE FROM THE SEGMENT START POINT TO THE POINT OF 
INTERSECTION OF EACH CROSSED OBSTACLE REFERENCE LINE. THESE 
DISTANCES ARE THEN SORTED FROM SMALLEST TO LARGEST AND CROSSINGS 

ARE UPDATED IN THE LINK LIST ACCORDINGLY. 


INPUT: HEAD() AND ALL COORDINATES REQUIRED TO DETERMINE THE 
DISTANCES FROM INITIAL VERTEX OF PATH SEGMENT TO EACH OF 
THE CROSSED REFERENCE LINES 


OUTPUT: HEAD(), NEXT() AND PRED() ARRAYS CONTAINING THE RAW STRING 
WHICH ACCURATELY REPRESENTS THE PATH ALONG THE CURRENT 
SEGMENT 


INTEGER SEG, LINE, HEAD(1000), NEXT( 1000), FSTSTR, HDINDX, 
+ START, FIRST( 1000), PRED(1000), F 
REAL*8 XINTER( 1000), YINTER( 1000), A1( 1000), B1(1000), 
+ A2( 1000), 82(1000), D2( 1000), DIST( 1000), X(1000), Y( 1000) 


DO 8 J = FSTSTR, LSTSTR 
LINE = ABS(HEAD(J)) 
CALL CASES2 (Al, B1, A2, B2, D2, SEG, LINE, XINTER, YINTER) 
DIST(J) = (XINTER(LINE)-X(SEG) )*¥*2+(YINTER(LINE)-Y (SEG) )**2 
DIST(J) = -DIST(J) 
NEXT(J) = J+1 
PRED(J+1) = J 

8 CONTINUE 


NEXT (LSTSTR) 


= 0 
PRED(LSTSTR+1) = 


LSTSTR 
CALL MERG2(DIST, NEXT, PRED, FSTSTR, F) 


NEXT(LOLD) = F 
PRED(F) = LOLD 
IF(SEG.EQ.1) START = LOLD 
LOLD = NEXT(LOLD) 
17 IF(NEXT(LOLD).NE.O) THEN 
LOLD = NEXT(LOLD) 
GO TO 17 
ENDIF 


RETURN 
END 
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SUBROUTINE CASES2(A1, B1, A2, B2, D2, SEG, LINE, XINTER, YINTER) 
EAR AEA EA EEE ER AEE RE EE EEE EE EEE 


DETERMINES COORDINATES FOR THE POINT OF INTERSECTION OF THE PATH 
SEGMENT AND EACH OF THE REFERENCE LINES IT CROSSES. AGAIN, IN ORDER 
TO MAINTAIN NUMERICAL STABILITY, CALCULATIONS ARE MADE USING THE 
MOST APPROPRIATE CLOSED FORM EQUATION. 


INPUT: COEFFICIENTS OF LINEAR EQUATIONS FOR CURRENT PATH SEGMENT 
AND ALL OBSTACLE REFERENCE LINES WHICH IT CROSSES 


OUTPUT: COORDINATES FOR POINTS OF INTERSECTION OF EACH OBSTACLE 
REFERENCE LINE WITH THE PATH SEGMENT 


IMPLICIT REAL*¥8 (A-H, O-Z) 


INTEGER SEG 

DIMENSION A1( 1000), B1( 1000), A2(1000), 82(1000), D2( 1000), 
+ XINTER( 1000), YINTER( 1000), A(2,2) 

A(1,1) = A1(LINE) 

A( 1,2) = BICLINE) 

A(2,1) = A2(SEG) 

A(2,2) = B2(SEG) 

BIGEST = 0.0 

DO 1L = 


V2 
DO 1K = 1,2 
TEST = DABS(A(K,L)) 

IF (TEST.GT.BIGEST) THEN 


BIGEST = TEST 


KBIG = K 
LBIG=L 
ENDIF 


1 CONTINUE 


IF (KBIG.EQ.1) THEN 

IF (LBIG.EQ.1) THEN 
YINTER(LINE)= -D2(SEG)/(B2(SEG)-B1(L INE )*A2(SEG)/A1(LINE)) 
XINTER(LINE)= -B1(LINE)*YINTER(LINE)/A1(LINE) 
RETURN 

ELSE 
XINTER(LINE)= -D2(SEG)/(A2(SEG)-A1(LINE)*#B2(SEG)/B1(LINE)) 
YINTER(LINE)= -A1(LINE)*XINTER(LINE)/B1(LINE) 
RETURN 

ENDIF 

ELSE 

IF (LBIG.EQ.1) THEN 

YINTER(LINE) = (D2(SEG)*A1(LINE)/A2(SEG) )/ 
+ (B1(LINE)-B2(SEG)*A1(L INE) /A2( SEG) ) 


68 


XINTER(LINE) = (-D2(SEG)-B2(SEG)*YINTER(LINE) )/A2(SEG) 


RETURN 

ELSE 
XINTER(LINE) = (D2(SEG)*B1(LINE)/B2(SEG) )/ 

+ (A1(LINE)-A2(SEG)*B1(LINE) /B2(SEG) ) 

YINTER(LINE) = (-D2(SEG)-XINTER(LINE)*A2(SEG) )/B2(SEG) 
RETURN 

ENDIF 

ENDIF 
END 


SUBROUTINE MERG2 (DIST, NEXT, PRED, FSTSTR, F) 
SRM SC RI IR HC CIE CIE SI CARE I TTI A RCE RC PCR FC CC C2 AR CK RC aC a a a ac aE a ic KC aR aC a aR a ae 


THIS SUBROUTINE SORTS A SUBSTRING OF ALL POSTIVE INTEGERS INTO 
INCREASING ORDER. 


INPUT: A DOUBLE LINK LIST CONSISTING OF DIST, NEXT, AND PRED ARRAYS 


OUTPUT: A DOUBLE LINK LIST WITH ALL ENTRIES PLACED IN THE ORDER 
IN WHICH THEIR RESPECT!]VE REFERENCE LINES WERE CROSSED 


IMPLICIT INTEGER(A-Z) 

REAL*8 DIST( 1000) 

LOGICAL DONE 

DIMENSION NEXT(1000), PRED( 1000) 


DONE = .FALSE. 
P = FSTSTR 


1F =P 
CALL SORT2 (F, PREDF, P, DIST, NEXT, DONE, PRED) 
IF(DONE) RETURN 
GOTO 1 


END 
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SUBROUTINE SORT2(F, PREDF, P, DIST, NEXT, DONE, PRED) 
RRR EEK EKER KR KERHK REAR KKH KERKAKE EA KK HK EK KKK KKEKKEREKKKKEKKKKEAKKAEKERE 


MARCHES DOWN A SUBSTRING OF POSITIVE INTEGERS AND ONLY SORTS IF 
ELEMENTS ARE NOT IN INCREASING ORDER. IF A NUMBER NEEDS TO BE 
PLACED HIGHER IN THE LIST SUBROUTINE ‘PUT’ 1S CALLED TO DO SO 


INPUT: POINTER F INTO DOUBLE LINK LIST ARRAYS DIST, NEXT, ANO PRED 
TO INDICATE THE BEGINNING OF A SUBSTING OF POSITIVE INTEGERS, 
AND PREDF 


OUTPUT: P IS A POINTER, DIST(P) IS THE LAST POSITIVE INTEGER IN THE 
SUBSTRING THAT IS BEGUN BY DIST(F) 


IMPLICIT INTEGER (A-Z) 

REAL*8 DIST( 1000) 

LOGICAL DONE 

DIMENSION NEXT( 1000), PRED( 1000) 
TAIL = F 

NTAIL = NEXT(TAIL) 


ls 


IF(NTAIL.EQ.0) THEN 
DONE = .TRUE. 
RETURN 

ENDIF 


IF(DIST(NTAIL).LE.DIST(TAIL)) THEN 
TAIL = NTAIL 
GOTO 1 
ENDIF 
CALL PUT2(F, TAIL, NTAIL, PREDF, DIST, NEXT, PRED) 
GOTO 1 


END 
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SUBROUTINE PUT2 (F, TAIL, NTAIL, PREDF, DIST, NEXT, PRED) 
C KEKKKKKKKKKKKKKKKKKKKKKKKKKKKEKKKKKKKKKKKKKKKKKKKKAKKKKAKKKEKKKKKKKKK ES 


C REARRANGES POINTERS TO PLACE POSITIVE INTEGERS IN INCREASING ORDER. 


INPUT: F START OF POSITIVE SUBSTRING 
TAIL - END OF CURRENTLY SORTED PORTION OF SUBSTRING 


NTAIL = NEXT(TAIL) - POINTER TO THE SUCCESSOR OF TAIL IN THE 
ITEM TO BE INCORPORATED INTO THE SORTED PORTION OF 
THE LIST 

PREDF - THE PREDECESSOR OF F IN LINKED LIST 


OUTPUT: SOME POINTERS IN NEXT AND PRED ARE CHANGED TO PUT DIST(NTAIL) 
IN ITS PROPER PLACE IN THE SORTED PORTION ON THE LIST 


OO AQOAN29N090N009 


IMPLICIT INTEGER (A-Z) 
REAL*8 D1IST( 1000) 
DIMENSION NEXT(1000), PRED( 1000) 


IF (DISTCNTAIL).GE.DIST(F)) THEN 
TEMP = NEXT(NTAIL) 
NEXT(NTAIL) = F 
PRED(F) = NTAIL 
NEXT(TAIL) = TEMP 
PRED(TEMP) = TAIL 
F = NTAIL 
RETURN 

ENDIF 

CALL WALK2(F, TAIL, NTAIL, DIST, NEXT, PRED) 

RETURN 

END 


viel 
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SUBROUTINE WALK2(F, TAIL, NTAIL, DIST, NEXT, PRED) 
SE RI IE ER SE EK RE IR A ICE I EE A IK EE SK FC I CK FR KE RAK ARIK KR HR 


THIS SUBROUTINE WALKS DOWN THE LINKED LIST AND PLACES THE DIST(NTAIL) 
IN THE CORRECT POSITION IN THE DIST ARRAY. 


INPUT: F, TAIL, NTAIL ARE USED AS ABOVE 


OUTPUT: ALTERS POINTERS IN NEXT AND PRED ARRAYS TO PLACE DIST(NTAIL) 
AFTER DIST(F) AND BEFORE DIST(TAIL) IN DOUBLE LINK LIST 


IMPLICIT INTEGER (A-Z) 
REAL*¥8 DIST( 1000), HNTAIL 
DIMENSION NEXT(1000), PRED( 1000) 


| =F 
NEXT! = NEXT(1) 
HNTAIL = DIST(NTAIL) 


os 


IFCHNTAIL .GE.DIST(NEXTI)) THEN 
NEXT(I) = NTAIL 


PRED(NTAIL) = | 
NNTAIL = NEXT(NTAIL) 
NEXT(NTAIL) = NEXTI 
PRED(NEXTI) = NTAIL 
NEXT(TAIL) = NNTAIL 
PRED(NNTAIL) = TAIL 
RETURN 

ENDIF 

| = NEXTI 

NEXTI = NEXT(1) 

GOTO 1 


END 


dz 


SUBROUTINE COUNTR(START, NEXT, HEAD, NELEMS) 
CRI CH RCE HE EE RE EE EE RR IE TE RE I EI IE IE RE RE I RE I I OK EK 


C THIS SUBROUTINE COUNTS THE NUMBER OF ELEMENTS IN A GIVEN STRING OF 
C CHARACTERS 


C INPUT: HEAD() AND NEXT() ARRAYS FOR THE STRING OF CHARACTERS 
C OUTPUT: NUMBER OF ELEMENTS IN THE STRING 
INTEGER PTR, NEXT(1000), HEAD(1000), START 


PTR = NEXT(START) 
NELEMS = 0 


20 IF (PTR.NE.O) THEN 
NELEMS = NELEMS + 1 
PTR = NEXT(PTR) 
GoriC, 20 
ENDIF 


RETURN 
END 


iS 
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SUBROUTINE ALG2(NOBS, NELEMS, HEAD, NEXT, X, Y, BX, BY, FR, NR, 
+ PR) 
HRA AAR ARATE ER EEE RATE AR KAR RAK ARERR RARER AA ART RTE 


GENERATES THE WELL-KNOW FUNDAMENTAL GROUP REPRESENTATION OF AN 
EQUIVALENCE CLASS. 


INPUT: RAW OR CANONICAL STRING, WITH # OBSTACLES, # ELEMENTS IN 
STRING, COORDINATES OF A AND OBSTACLES (BK) 


OUTPUT: 1) IF INPUT 1S RAW STRING, THEN POSSIBLY UNREDUCED 
FUNDAMENTAL GROUP WILL RESULT (1F CANCELLATION OF LIKE 
POS!T!VE NUMBERS COULD HAVE OCCURRED IN THE RAW STRING). 
2) IF INPUT 1S CANONICAL STRING, THEN THE RESULTING 
FUNDAMENTAL GROUP WILL BE IN REDUCED FORM. 


REAL*8 XA, YA, XB, YB, X(1000), Y( 1000), BX(1000), BY( 1000) 
INTEGER NOBS, NELEMS, HEAD( 1000), NEXT( 1000), PRED( 1000), 
+ FR( 1000), NR( 1000), PR( 1000), $( 1000), FUNDGP( 1000), 
+ START 

LOGICAL RIGHT( 1000), RITE 


LENGTH = 1 
START = NEXT(1) 


X(1) 
i) 


DO 1K = 1, NOBS 
XB = BX(K) 
YB = BY(K) 
1 RIGHT(K) = RITE(XA, YA, XB, YB) 


XA 
YA 


NELEM = NELEMS + 1 
DO 10 M = 2, NELEM 
S(M) = HEAD(START) 
J = 1ABS(S(M)) 
RIGHT(J) = .NOT. RIGHT(J) 
IF(S(M).GT.0) THEN 
LENGTH = LENGTH + 1 
IF(RIGHT(J)) THEN 


FUNDGP(LENGTH) = J 
EESE 
FUNDGP (LENGTH) = -J 
ENDIF 
ENDIF 


START = NEXT(START) 
IF(START.EQ.0) GOTO 11 
10 CONTINUE 


74 


HeERCI), = 0 
NR(1) = 2 
DO 2 | = 2, LENGTH 
FR(I) = FUNDGP(1) 
NR(I) = 1 + 1 
PR(I) = 1 - 1 
2 CONTINUE 


NR(LENGTH) = 0 


RETURN 
END 


LOGICAL FUNCTION RITE(XA, YA, XB, YB) 
HEAR A AER ACA AE EEE AE AE I AE AC AC AE AC AE AEE A AE AE AC AE A A A AE A IE AC A CO CAR A A A A A A A A AC AC A OR AO AR OR A a 


THIS FUNCTION DETERMINES WHICH SIDE OF A GIVEN DIRECTED LINE ANY 
POINT LIES. 


INPUT: TWO POINTS THAT DETERMINE THE LINE 


OUTPUT: LOGICAL VARIABLE THAT 1S TRUE IF A POINT LIES TO THE RIGHT 
AND FALSE IF A POINT LIES TO THE LEFT 


REAL*8 XB, YB, XA, YA, SIGNA 
RITE = .TRUE. 
SIGNA = -((YB*¥XA)-(XB*YA) ) 
IF(SIGNA.LT.O)THEN 
RITE = .FALSE. 
RETURN 
ENDIF 
IF(SIGNA.EQ.0) THEN 
PRINT¥, ‘THE POINT A LIES ON THE LINE LK PROGRAM STOPS' 
STOP 
ENDIF 


RETURN 
END 
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SUBROUTINE CANALG2(HEAD, NEXT, PRED) 
SE AACA IRR EE EE A HE EE EEE CE AE EE EE EK EE ER A OR AE OK AE A A EE EK OE OR EK AK OK 


THIS SUBROUTINE TAKES AN UNREDUCED FUNDAMENTAL GROUP REPRESENTATION 
OF A GIVEN CLASS AND CANCELS A GENERATOR IF IT IS ADJACENT TO ITS 
INVERSE. 


INPUT: HEAD, NEXT, AND PRED ARRAYS 


OUTPUT: HEAD, NEXT, AND PRED ARRAYS WITH NEXT AND PRED REARRANGED 
TO SKIP AROUND CANCELLED ELEMENTS 


IMPLICIT INTEGER (A-Z) 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


START = 1 

PTR1 = START 
PTR2 = START 

PTR3 = NEXT(PTR2) 


10 IF (PTR3.NE.0) THEN 
IF (HEAD(PTR2).EQ.-(HEAD(PTR3))) THEN 
NEXT(PTR1) = NEXT(PTR3) 
PRED(NEXT(PTR3)) = PTR1 
\F(NEXT(PTR1).EQ.0) RETURN 


PTR2 = PTR1 
PTR1 = PRED(PTR1) 
PTR3 = NEXT(PTR2) 
ELSE 
PTR1 = PTR2 
PTR2 = PTR3 
PTR3 = NEXT(PTR2) 
ENDIF 
GO TO 10 
ENDIF 
RETURN 
END 
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SUBROUTINE ALG1(HEAD, NEXT, PRED, NELEMS) 
AOR RCI IR IR RC AR CR RR AE ROR IR A RO OR A RAR I RRC I I RR A A ROR A AOR ROR RO OK I 


THIS SUBROUTINE TAKES A GIVEN RAW STRING OF CHARACTERS REPRESENTING 
A PATH AND PRODUCES THE CANONICAL FORM OF THAT STRING. THIS IS DONE 
BY FIRST ORDERING ALL OF THE ALPHA SUBSTRINGS FROM SMALLEST TO 
LARGEST IN ABSOLUTE VALUE. NEXT, CANCELLATION IS PERFORMED TO 
ELIMINATE ALL LIKE PAIRS OF ADJACENT ELEMENTS FROM THE STRING. 


INPUT: RAW STRING IN FORM OF DOUBLE LINKED LIST WITH HEAD(), NEXT() 
AND PRED() ARRAYS 


OUTPUT: CANONICAL FORM OF THE RAW STRING 


IMPLICIT INTEGER (A-Z) 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


START = 1 


CALL MERGI(HEAD, NEXT, PRED) 
CALL CANCEL(START, HEAD, NEXT, PRED, NELEMS) 


RETURN 
END 


of 


OO 
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SUBROUTINE MERG1 (HEAD, NEXT, PRED) 
ERK KKK KK HK HK HK EK EK KE HE RE RK EK ER EK EK EK EK EE KK AE OK KE OK I OK EK KK IK KKK 


THIS SUBROUTINE DOES THE INITIAL SORT OF THE ALPHA SUBSTRINGS IN THE 
RAWSTRING (INCREASING IN ABSULUTE VALUE). 


INPUT: HEAD, NEXT, PRED ARRAYS REPRESENTING A DOUBLE LINKED LIST 
OF THE RAWSTRING WITH THE ALPHA SUBSTRINGS UNORDERED 


OUTPUT: POINTERS STORED IN THE ARRAYS NEXT AND PRED ARE ALTERED SO 
THAT EACH SUBSTRING OF THE STORED LIST WHICH CONSISTS 
ENTIRELY OF NEGATIVE INTEGERS IS SORTED INTO NON-INCREASING 
ORDER, WHILE SUBSTRINGS OF POSITIVE INTEGERS ARE LEFT 
UNALTERED. 

IMPLICIT INTEGER(A-Z) 
LOGICAL DONE 
DIMENSION HEAD( 1000), NEXT(1000), PRED( 1000) 


DONE = .FALSE. 
pa 


1 CALL FRONT (P, F, PREDF, HEAD, NEXT, DONE) 
IF (DONE) RETURN 
CALL SORT (F, PREDF, P, HEAD, NEXT, DONE, PRED) 
|F (DONE) RETURN 
GOTO 1 


END 
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SUBROUTINE FRONT(P, F, PREDF, HEAD, NEXT, DONE) 
C FRERKKKKKKKKKKKKAKKKKKAKKAKKAKAKKAKKAKKAK AKA KAKAKKKAKAKKAKKKKKKKAKAKAKKKAKKKK KK KE 


C FINDS THE BEGINNING OF NEGATIVE INTEGER STRINGS (ALPHA STRING) 
C INPUT: POINTER P INTO LINKED LIST 
OUTPUT: F IS A POINTER INTO THE HEAD ARRAY. POINTS TO FIRST NEGATIVE 


ENTRY WHICH OCOURS STRICTLY AFTER HEAD(P). PREDF IS POINTER 
SUCH THAT HEAD(PREDF) 1S THE PREDESSOR OF HEAD(F). 


OOO 


IMPLICIT INTEGER(A-Z) 
LOGICAL DONE 
DIMENSION HEAD( 1000), NEXT( 1000) 


Ea = uP 
1 PREDF = F 

F = NEXT(F) 

IF (F.EQ.0) THEN 
DONE = .TRUE. 
RETURN 

ENDIF 


IF (HEAD(F).LT.O) RETURN 
GOTO 1 


END 
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SUBROUTINE SORT(F, PREDF, P, HEAD, NEXT, DONE, PRED) 


TEER OR IIE IC IE 2 IE 2 IC aE AIC I ICI 2 a IC a AS AE IC IE RE I AR RC OE IE RC IC aR A I RE OR OE RR RR OR 


MARCHES DOWN A SUBSTRING OF NEGATIVE INTEGERS AND ONLY SORTS IF 
ELEMENTS ARE IN INCREASING ORDER. IF A NUMBER NEEDS TO BE 
PLACED HIGHER IN THE LIST SUBROUTINE ‘PUT’ IS CALLED TO DO SO 


INPUT: POINTER F INTO DOUBLE LINK LIST ARRAYS HEAD, NEXT, AND PRED 


TO INDICATE THE BEGINNING OF A SUBSTRING OF NEGATIVE 
INTEGERS; AND PREOF 


OUTPUT: P IS A POINTER, HEAD(P) IS THE LAST NEGATIVE INTEGER IN THE 


a 


SUBSTRING THAT IS BEGUN BY HEAD(F) 


IMPLICIT INTEGER (A-Z) 
LOGICAL DONE 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


TAIL = F 
NTAIL = NEXT(TAIL) 


IF(NTAIL.EQ.0) THEN 
DONE = . TRUE. 
RETURN 

ENDIF 


IF (HEAD(NTAIL).GT.O) THEN 
P = TAIL 
RETURN 
ENDIF 
IF(HEAD(NTAIL).LE.HEAD(TAIL)) THEN 
TAIL = NTAIL 
GOTO 1 
ENDIF 
CALL PUT(F, TAIL, NTAIL, PREDF, HEAD, NEXT, PRED) 
GOTO 1 


END 
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YOAOANI909N90 9090 


OO 


SUBROUTINE PUT (F, TAIL, NTAIL, PREDF, HEAD, NEXT, PRED) 
FRR ROR FE FR OI OI TRIO I I IE ICC IIE AE RRC RCE RCC RCE RCE RCE AIC AE ICC RCE RCC ICC IK RC ICC i I RE IE OE 


SUBROUTINE THAT REARRANGES POINTERS TO PLACE A NEGATIVE INTEGER IN 
NON- INCREASING ORDER. 


INPUT: F ~ START OF NEGATIVE SUBSTRING 
TAIL END OF CURRENTLY SORTED PORTION OF SUBSTRING 


NTAIL = NEXT(TAIL) - POINTER TO THE SUCCESSOR OF TAIL IN THE 
ITEM TO BE INCORPORATED INTO THE SORTED PORTION OF 
tek oil. 

PREDF - THE PREDECESSOR OF F IN LINKED LIST 


OUTPUT: SOME POINTERS !N NEXT AND PRED ARE CHANGED TO PUT HEAD(NTAIL) 
IN !TS PROPER PLACE IN THE SORTED PORTION ON THE LIST 


IMPLICIT INTEGER (A-Z) 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


IF (HEAD(NTAIL).GE.HEAD(F)) THEN 
NEXT(PREDF) = NTAIL 
PRED(NTAIL) = PREDF 
TEMP = NEXT(NTAIL) 
NEXT(NTAIL) = F 
PRED(F) = NTAIL 


NEXT(TAIL) = TEMP 
PRED( TEMP) = TAIL 
F = NTAIL 
RETURN 

ENDIF 


CALL WALK(F, TAIL, NTAIL, HEAD, NEXT, PRED) 


RETURN 
END 
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SUBROUTINE WALK(F, TAIL, NTAIL, HEAD, NEXT, PRED) 
CREE IC EEE EE EE HE EEE EE EE EEE EE EEE RE EEE ER RR OR IO OOK ak 


C THIS SUBROUITNE WALKS DOWN THE LINKED LIST AND PLACES THE HEAD(NTAIL) 
C IN THE CORRECT POSITION IN THE HEAD ARRAY. (DECREASING ORDER) 


C INPUT: F, TAIL, NTAIL ARE USED AS ABOVE 


C OUTPUT: ALTERS POINTERS IN NEXT AND PRED ARRAYS TO PLACE HEAD(NTAIL ) 
C AFTER HEAD(F) AND BEFORE HEAD(TAIL) !N DOUBLE LINK LIST 


IMPLICIT INTEGER (A-Z) 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


| =F 
NEXT! = NEXT(1) 
HNTAIL = HEAD(NTAIL) 


1 IFCHNTAIL.GE.HEAD(NEXTI)) THEN 
NEXT( I) = NTAIL 
PRED(NTAIL) = | 
NNTAIL = NEXT(NTAIL) 
NEXT(NTAIL) = NEXTI 
PRED(NEXTI) = NTAIL 


NEXT(TAIL) = NNTAIL 
PRED(NNTAIL) = TAIL 
RETURN 

ENDIF 

| = NEXTI 

NEXT| = NEXT(1) 

GOTO 1 


END 
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SUBROUTINE CANCEL(START, HEAD, NEXT, PRED, NELEMS) 


HM RK AHHH HH AE RC IC AE iC ICR RIC A A HG ACK ARC AC fC IC IC RG RC AIC IC IC ICRC AIC IC IC RC RC IR A AIC OR A AO OR RO OK 


GIVEN THE SORTED RAW STRING, THIS SUBROUTINE REDUCES THE STRING TO 
CANONICAL FORM BY ADJUSTING THE LINK LIST TO SKIP ANY PAIRS OF 
ADJACENT LIKE ELEMENTS IN THE STRING. ONCE TWO ELEMENTS ARE 
ELIMINATED FROM THE STRING, THE SUBSTRINGS WHICH WERE ON THEIR LEFT 
AND RIGHT CONCATENATE TO FORM A NEW STRING AND HENCE A NEW PAIR OF 
ADJACENT ELEMENTS WHICH MUST ALSO BE CHECKED FOR EQUALITY. IN THE 
EVENT A BETA SUBSTRING BECOMES ANNHILATED FROM THE STRING, THE TWO 
ADJACENT ALPHA STRINGS CONCATENATE AND ARE AGAIN SORTED FROM SMALLEST 
TO LARGEST IN ABSOLUTE VALUE AS A SINGLE SUBSTRING. 


INPUT: SORTED RAW STRING 
OUTPUT: CANONICAL FORM OF THE RAW STRING 


IMPLICIT INTEGER (A-Z) 
DIMENSION HEAD( 1000), NEXT( 1000), PRED( 1000) 


START = 1 
START1 = START 
START2 = START 


BETA = START 
PTR1 = START1 
PTR2 = START1 

50 PTR3 = NEXT(PTR2) 


IF (HEAD(PTR3).LT.0) THEN 
START1 = PTR3 
START2 = START1 

BESE 
BETA = PTR3 

ENDIF 


10 IF (PTR3.NE.0) THEN 
CALL CHEKER (HEAD, START1, START2, PTR1, PTR2, PTR3, BETA) 
IF CHEAD(PTR2) .EQ.HEAD(PTR3)) THEN 
NEXT(PTR1) = NEXT(PTR3) 
PRED(NEXT(PTR3)) = PTR1 
PTR2 = NEXT(PTR3) 
IF (PTR2.EQ.0) GO TO 60 
CHECK 1 = HEAD(PTR2) * HEAD(PTR3) 
CHECK2 = HEAD(PTR1) * HEAD(PTR2) 
IF (CHECK1.LT.0.AND.CHECK2.GT.0.AND.HEAD(PTR2).LT.O)THEN 
PTR2 = PTR1 
PTR1 PRED(PTR1) 
PTR3 = NEXT(PTR2) 
START1 = START2 
6 IF (START1.NE.1.AND.HEAD(PRED(START1)).LT.0O)THEN 
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START1 = PRED(START1) 
GO TO 6 
ENDIF 
START2 = PTR3 
CALL MERGE (HEAD, NEXT, START1, START2, PRED) 
PTR1 = PRED(START1) 
PTR2 = PTR1 
GO TO 50 
ELSEIF (CHECK1.LT.0.AND.CHECK2.EQ.0.AND.HEAD(PTR2).LT.0) 
+ THEN 
PTR2 = PTR1 
PTR3 = NEXT(PTR1) 
START1 = PTR3 
START2 = START1 
GO TO 10 
ELSEIF (CHECK1.GT.0.AND.CHECK2.GT.O.AND.HEAD(PTR2).GT.0) 
+ THEN 
IF (BETA.EQ.2) THEN 
BETA = 1 
ELSE (BETA.EQ 41) alneEN 
GO TO 9 
ELSE 
IF(START1.NE.1) BETA = PRED(START1) 
11 IF (HEAD(PREO(BETA)).GT.0) THEN 
BETA = PRED(BETA) 
Go TO a 
ENDIF 
ENDIF 
3 PIR 4 


Wott 


PREO( BETA) 
PTR2 = BETA 
PTR3 = NEXT(PTR2) 
GO TO 10 
ELSEIF (CHECK1.LT.0.AND.CHECK2.GT.0.AND.HEAD(PTR2) .GT.0O) 
te THEN 
PTR2 
PTR 1 


PTR 
PRED(PTR1) 


IF (HEAD(PTR1).LT.0) THEN 
START1 = PTR1 
> IF (HEAD(PRED(START1)).LT.O)THEN 
START1 = PRED(START1) 
GO Ter 
ENDIF 
ENDIF 
PTR3 = NEXT(PTR2) 
GO TO 10 
ELSEIF (CHECK1.LT.0.AND.CHECK2.LT.O.AND.HEAD(PTR2).LT.O) 
+ THEN 
IF (START1.EQ.1) THEN 
START1 = PTR2 
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START2 = PTR2 
PTR3 = NEXT(PTR2) 
ELSE 
START1 = STARTZ 
PTR3 = NEXT(PTR2) 
ENDIF 
GO TO 10 
ELSE 
START1 = START2 
PTR3 = NEXT(PTR2) 


ENDIF 
ELSE 
PTR1 = PTR2 
PTR2 = PTR3 
PTR3 = NEXT(PTR2) 
ENDIF 
GO TO 10 
ENDIF 


60 CALL COUNTR (START, NEXT, HEAD, NELEMS) 


RETURN 
END 
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SUBROUTINE CHEKER(HEAD, START1, START2, PTR1, PTR2, PTR3, BETA) 
AER RE RE EEE HE AE EO RE OE EE IEE EE IEE EHEC IC EEA CAR IC IE IK IK CC a a a aK 2K aK a aK a a 2K a aK ok 2 


AS POINTERS MOVE ALONG THE CHARACTER STRING DURING CANCELLATION THIS 
SUBROUTINE CHECKS TO DETERMINE WHETHER OR NOT THE END OF ONE 
SUBSTRING IS REACHED AND A NEW ONE BEGINS. DEPENDING ON THE OUTCOME 
OF THIS CHECK, THE POINTERS USED TO IDENTIFY THE BEGINNING OF THE 
TWO LATEST ALPHA STRINGS AND THE LATEST BETA STRING ARE UPDATED 


INPUT: HEAD() ARRAY AND PRESENT POINTER LOCATIONS FROM ‘CANCEL ' 
SUBROUT | NE 


OUTPUT: MODIFIED INDICES FOR LOCATION OF ALPHA AND BETA STRINGS 
INTEGER HEAD( 1000), START1, START2, PTR1, PTR2, PTR3, CHECK, BETA 
CHECK = HEAD(PTR3) * HEAD(PTR2) 


1[F (CHECK.GT.0) THEN 
1F (HEAD(PTR2).LT.O.AND.HEAD(PTR1).GT.0) START2 
[F (HEAD(PTR2).GT.O0.AND.HEAD(PTR1).LT.0O) BETA = 
ELSEIF (CHECK.LE.0O) THEN 
IF (HEAD(PTR2).GT.O.AND.HEAD(PTR1).LT.0O) THEN 
BETA = PTR2 
START1 = START2 
START2 = PTR3 
ELSEIF (HEAD(PTR2).GT.0) THEN 
If CSMART 1.NE. 1) ° THEN 
START1 = START2 
START2 PTR3 
ELSE 
START 1 
START 2 
ENDIF 
ELSEIF (HEAD(PTR2).LT.O.AND.HEAD(PTR1).GT.0) THEN 
BETA = PTR3 
START1 = START2 
START2 = PTR2 
ELSEIF (HEAD(PTR2).LT.O.AND.HEAD(PTR1).LT.0) THEN 
BETA = PTR3 
ENDIF 
ENDIF 


= PTR2 
PTR2 


PTR3 
START 1 


RETURN 
END 
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te: 


SUBROUTINE MERGE(HEAD, NEXT, START1, START2, PRED) 
HAHAH AHR HAE HA RA A HC HE AE AE HR AE HE HO RAC A AE AC AG ACE AC RCA A AR AC ACE AOR AR AOE AIK AA A RO oe 


THIS SUBROUTINE MERGES TWO ORDERED ALPHA STRINGS. THE RESULTING 
SUBSTRING WILL BE INCREASING IN ABSOLUTE VALUE. 


INPUT: HEAD(), NEXT() AND PRED() ARRAYS ALONG WITH INDICES START1 
AND START2 INDICATING BEGINNING OF THE ALPHA STRINGS TO BE 
MERGED 


OUTPUT: MODIFIED STRING CONTAINING THE MERGED ALPHA STRINGS WHICH 
WERE A RESULT OF ANNHILATION OF A BETA STRING 
INTEGER HEAD( 1000), NEXT( 1000), START1, START2, P1, P2, TAIL, 
+ PRED( 1000 ) 
LOGICAL GO 


IF (START1.EQ.1) GO TO 2 
IF (HEAD(PRED(START1)).EQ.0) START1 = NEXT(PRED(START 1) ) 


|F (NEXT(PRED(START1)).NE.START1)THEN 
START1 = NEXT(PRED(START1) ) 
i IF (START1.EQ.1) GO TO 2 
|F (HEAD(START1).GT.0) THEN 
START1 = NEXT(START1) 
GO TO 7 
ENDIF 
8 |F (HEAD(PRED(START1)).LT.0O)THEN 
START 1 = PRED(START1) 
GO TO 8 
ENDIF 
ENDIF 


2 CALL SET(HEAD, NEXT, START1, START2, P1, P2, TAIL, PRED) 
1 IF(GO(P1, P2, START2, HEAD, NEXT)) THEN 
CALL MERGER(HEAD, NEXT, P1, P2, TAIL, PRED) 
GOTO 1 
ENDIF 
CALL FXTAIL(HEAD, NEXT, P1, P2, TAIL, START2, PRED) 


RETURN 
END 
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SUBROUTINE SET(HEAD, NEXT, START1, START2, P1, P2, TAIL, PRED) 
HEA RAE RAE HE REE A HE RAE RC CE RAE A AE ICA A IC ICA A RC I RA ICR ICRA CIA AE IORI A A ICA A A A A A a 


SETS THE STARTING POINTER OF THE NEW STRING AND THE TWO POINTERS 
OF THE STRINGS TO BE MERGED 


INPUT: POINTERS INTO TWO ALPHA STRINGS 
OUTPUT: POINTER INTO BEGINNING OF A NEW SORTED STRING 


INTEGER HEAD( 1000), NEXT( 1000), START1, START2, P1, P2,TAIL, 
+ PRED( 1000) 


P1 = START1 

P2 = START2 

IF (HEAD(P1).GT.HEAD(P2)) THEN 
TAIL = P1 
START1 = P1 
P1 = NEXT(P1) 

ELSE 
TAIL = P2 
START1 = P2 
NEXT(PRED(P1)) = P2 
PRED(P2) = PRED(P1) 
P2 = NEXT(P2) 

ENDIF 


RETURN 
END 


88 


LOGICAL FUNCTION GO(P1, P2, START2, HEAD, NEXT) 
CR RE EK HK KK EE HE EE HE EE EE EE EE EE EE AE A A OK CAE AE OK OE AE CK OK AE AE EK OK 


C LOGICAL FUNCTION TO DETERMINE WHEN THE END OF EITHER LIST IS REACHED 
C INPUT: POINTERS INTO TWO ALPHA SUBSTRINGS 


C OUTPUT: LOGICAL VARIABLE WHICH IS TRUE IF EITHER POINTER IS AT THE 
C END OF A STRING, AND FALSE OTHERWISE 


INTEGER HEAD( 1000), NEXT( 1000), P1, P2, START2 
GO = .TRUE. 


IF (P1.EQ.START2) THEN 
GO = .FALSE. 
RETURN 

ENDIF 


IF (NEXT (P2) .EQ.0)THEN 
GO = .FALSE. 
RETURN 

ENDIF 


IF (HEAD(P2).GT.0)THEN 
GO = .FALSE. 
RETURN 

ENDIF 


RETURN 
END 
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SUBROUTINE MERGER(HEAD, NEXT, P1, P2, TAIL, PRED) 
CRRA RE EE ER EK EK HE OK EE EO EA EE EE OK EE EO OE EO Eo 2K oo KC a OE IK OK 2 2K 2 2K OK 


THIS SUBROUTINE COMPARES TWO ELEMENTS - ONE IN EACH ALPHA SUBSTRING 
AND PLACES THE ONE THAT IS THE SMALLEST ON THE LIST REPRESENTING THE 
NEW SORTED LIST 


INPUT: HEAD, NEXT, PRED, AND POINTERS INTO THE HEAD ARRAY 


OUTPUT: HEAD, NEXT, PRED, AND POINTERS FOR THE NEXT TWO ELEMENTS 
TO BE CHECKED 


OL. “On 0 Ore 


INTEGER HEAD( 1000), NEXT(1000), P1, P2, TAIL , PRED( 1000) 


IF (HEAD(P1).GT.HEAD(P2)) THEN 

CALL TACON1(NEXT, P1, TAIL, PRED) 
EESE 

CALL TACON2(NEXT, P2, TAIL, PRED) 
ENDIF 


RETURN 
END 


SUBROUTINE TACON1 (NEXT, P1, TAIL, PRED) 
CMA HE EE KE EE KE AC ICC OK OE OE OE 2 OE OE OR A 2 OR OO oR KK He Ke eK ee Ee KC KC ee 2 2 2 2 2K 2K 2K 2K KK 


C ADDS THE ELEMENT TO WHICH P1 POINTS TO THE TAIL OF THE NEW STRING 


C INPUT: NEXT, PRED, Pl, TAlE 


© 


OUTPUT: NEXT, PRED, P1, TAIL 
INTEGER NEXT( 1000), P1, TAIL, PRED( 1000) 


NEXT(TAIL) = P1 
PRED(P1) = TAIL 
TAIL = P1 

P1 = NEXT(P1) 


RETURN 
END 


SUBROUTINE TACON2 (NEXT, P2, TAIL, PRED) 
C KERR KAKA HK KKAK KAKA HK KK AK AKA KA KAA HK AKA ARK AHH A KAKA HA KHER HHH EK H 


C ADDS THE ELEMENT TO WHICH P2 POINTS TO THE TAIL OF THE NEW STRING 


© INPOT: NEAT, PRED, PZ, TAIL 
C OUTPUT: NEXT, PRED, P2, TAIL 


INTEGER NEXT( 1000), P2, TAIL, PRED( 1000) 


NEXT(TAIL) = P2 
PRED(P2) = TAIL 
TAIL = P2 

P2 = NEXT(P2) 
RETURN 

END 
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SUBROUTINE FXTAIL(HEAD, NEXT, P1, P2, TAIL, START2, PRED) 
KHKKKKKKKKKKKKKKK KKK KKKKKKKKAKKKKKKKKKKKKKK KKK KK KKKKKKKKKKKKKKKKKKKKKKKSE 


ATTACHES THE END OF THE LONGER SORTED STRING TO THE END OF THE MERGED 


STRING 


INPUT: 


P1, P2, HEAD, NEXT, PRED, TAIL 


OUTPUT: P1, P2, HEAD, NEXT, PRED, TAIL, START2 


INTEGER HEAD(1000), NEXT(1000), P1, P2, TAIL, START2, P11, 


a 


PRED( 1000) 


IF(P1.EQ.START2) THEN 


CALL TACON2(NEXT, P2, TAIL, PRED) 
RETURN 


ENDIF 


1F (NEXT (P2).EQ.0)THEN 


2 


CALL TACON1(NEXT, P1, TAIL, PRED) 
P11 = TAIL 
IF (P1.NE.START2) THEN 
P11 = P1 
P1 = NEXT(P1) 
GOTO 2 
ENDIF 
NEXT(P11) = P2 
PRED(P2) = P11 
RETURN 


ENDIF 


IF (HEAD(P2) .GT.0)THEN 


3 


CALL TACON1(NEXT, P1, TAIL, PRED) 
ah) = UAL 
IF(P1.NE.START2) THEN 
P11 = PI 
P1 = NEXT(P1) 
GOTO 3 
ENDIF 
NEXT(P11) = P2 
PRED(P2) = P11 
RETURN 


ENDIF 


RETURN 


END 
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SUBROUTINE CTREX(FL, NL, PL, FR, NR, PR, BSEED, PSEED) 
CRRA IE KK AEE OK AK EK EE EK IK OK AE AE IC OK EK EK KC EAE AE OE EK AE AE IE AE EE OK AE EE 2K OK EE IE KC OE OK OK OK 


C THIS SUBROUTINE CHECKS F(C(R(P))) AGAINST K(FC(R(P))) FOR EQUALITY. 
C INPUT: FUNDAMENTAL GROUP REPRESENTATIONS FOR BOTH ALGORITHMS 


C OUTPUT: MESSAGE INDICATING THE OCCURRENCE OF UNEQUAL FUNDAMENTAL 
C GROUPS 


INTEGER FL( 1000), NL( 1000), PL( 1000), FR( 1000), NR( 1000), PR( 1000) 
IF(NL(1).EQ.0.AND.NR(1).EQ.0) GO TO 11 


J 


NL(1) 
NR( 1) 


10 IF (1.EQ.0.AND.J.EQ.0) GO TO 11 
IF (FL(1).EQ.0.AND.FR(J).EQ.0) GO TO 11 


IF (FLCI1).EQ.FR(J)) THEN 

| = NLCI) 

J = NR(J) 

GOTO 10 

ELSE 

PRINT*, STRINGS UNEQUAL ' 

PRINT¥,° — 

IF(NL(1).EQ.0)THEN 
PRINT*, ‘THE GROUP IS EMPTY AFTER ALG1 SEQUENCE’ 

ELSE 
PRINT*, THE FUNDAMENTAL GROUP AFTER ALG1 IS: 
CALL PRINTS(FL,NL) 

ENDIF 

IF(NR(1).EQ.0) THEN 
PRINT*, ‘THE GROUP IS EMPTY AFTER ALG2 SEQUENCE’ 

ELSE 
PRINT*, ‘THE FUNDAMENTAL GROUP AFTER ALG2 IS:' 
CALL PRINTS(FR,NR) 

ENDIF 

PRINT¥,° ¢ 

PRINT*, ‘BSEED AND PSEED ARE‘ ,BSEED,PSEED 


ENDIF 


11 RETURN 
END 
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SUBROUTINE PRINTS(HEAD, NEXT) 
HEA AA HAE AR IO ACA CIE AE AE ICE FC IE I CE A EA A A AR AR A IR AR RR RR OR OR AR RR a OK a ok oR aR ok ak 


THIS SUBROUTINE CAN BE USED TO PRINT OUT ANY STRING THAT !S STORED IN 
LINK LIST FORM 


INPUT: LINK LIST ARRAYS 
OUTPUT: HORIZONTAL STRING OF THE ELEMENTS ACCORDING TO THE NEXT() 
INTEGER HEAD( 1000), NEXT( 1000), STRING( 1000) 


NSTART = NEXT(1) 
| = Q 


98 IF(NSTART.NE.O) THEN 
pis ol ees i 
STRING(!) = HEAD(NSTART) 
NSTART = NEXT(NSTART) 
GOTO 98 
ENDIF 


PRINT 11152 (StRINGRI, 2 = 91,4) 
111 FORMAT(' *,2014) 
PRINT*,° 


RETURN 
END 
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ALL SUBROUTINES BELOW THIS POINT WERE USED AS DEBUG TOOLS. THEY 
GENERATE CRUDE GRAPHS WHICH PLOT THE BOARDS AND THE PATHS. 
THE VALUES AQUIRED BY THE RAWSTRING , MERGE, AND MANY OTHER 


SUBROUTINES WERE CHECKED USING THESE ROUTINES. 
A I OR A I RE EE OE OK IC I OK 2 OK SIE 2 IE a oC i 2 aC 2 2 iC 2 i oC I 2 ok ik Ci 2 ok i 2 a a ak a a ok 


OOAQ0N00 


SUBROUTINE GRAPH(BX, BY, X, Y, NOBS, NUMPTS) 
CR RR A RK AR EK EE A AE OK OK OK A OK OK KE KE OK OE OE EEK OE EE OK OE OK EE 2K OK KK KK 


REAL*8 X( 1000), Y(1000), Bx( 1000), BY( 1000) 
CHARACTER*¥1 MATR(53, 105) 


CALL GRAFPA(X, Y, NUMPTS, MATR) 
CALL GRAFOB(BX, BY, NOBS, MATR) 


pomin! = 1,53 
1 PRINT 111, (MATR(I,J), J=1, 105) 
111 FORMAT(' ',105A1) 


RETURN 
END 


vhs 


SUBROUTINE GRAFPA (X, Y, NUMPTS, MATR) 
CC KRKKKKKKKKKKAKKKA KRHA KAAK AKA A KAKA KKKKHA KKK KKAKH KA KKK KH KK AK KKK HK KKK KK KKK AK 


REAL¥*8 X( 1000), Y( 1000) 
CHARACTER*1 MATR(53, 105), CH 


DO i Beh sy 
10 M = 2,104 
10 MATR(L,M) = °° 
pO Tie) 105 
MATR(1,1) = '' 
1 MATR(53,1) = ‘*' 
DO Maa a= sams 
MATR(1,1) = 
4 MATR(1,105) = 
xa) 
YK= 1G) 


CALL COOR(XK, YK, IX, IY) 
CALL CHARPA(1, CH) 
MATRCIY,1X) = 


DO 2 K = 2,NUMPTS 
2 CALL FILL(K, X, Y, MATR) 


DO 3 K = 2,NUMPTS 
XK = X(K) 
YK = Y(K) 


CALI COOR TKK, XK ale) 
CALL CHARPA(K, CH) 
3 MATR(IY,1X) = 


RETURN 
END 


SUBROUTINE COOR(BX, BY, IX, IY) 
Co KKK KKKKKK KEK HKK KKK KKK KKK KA KKK KKK KKK HK KKK KKK KKK KKK KKK KKK KK KKKK HK KKKKKKKK 


REAL¥8 BX, BY 


IX = MESHP(BX) * 2 
\Y = 54 - MESHP(BY) 
RETURN 

END 
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INTEGER FUNCTION MESHP(X) 
CC FRR REAR AE AK RK RK OK AE OK OK RO EE EOE 2K A RO AEE AE OE A OR EE AE EOE EE OK I A IE 2 OK IK 2K 2 2K IK aK 2K 2K OK ok 


REAL*8 X 
MESH = C(Xtt. )#(51./72.)) + 2, 


RETURN 
END 


SUBROUTINE CHARPA(K, CH) 
CRA IER EK OR RE OK OK RR OK OR OK OE AE OK KC OK KC AE AE OK OK CA AE AC OE AE AE HE CA OE A AE AR HE AE OK A AE KOK OK AE OK 


CHARACTER* 1 CH 


IF(K.EQ.1) THEN 
CH = ‘A’ 
RETURN 

END IF 


IF(K.EQ.2) THEN 
Cag. Bi 
RETURN 

END IF 


IF(K.EQ.3) THEN 
CH= ‘C 
RETURN 

ENDIF 


(P(K.EQ.4) THEN 
Cie 
RETURN 

ENDIF 


IF(K.EQ.5) THEN 
CHl== E> 
RETURN 

ENDIF 


IF(K.EQ.6) THEN 
CH = ‘A’ 
RETURN 

ENDIF 


END 
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SUBROUTINE FILL(K, X, Y, MATR) 
C KKK KKKKKKKKKKKKKKKKK KKK KKKAK KKK KKKKHK KKK KKKKKKKKKKKKKKKKKKKAKKKKKAKKKKKKAE 


REAL*8 X( 1000), Y(1000), U, V, DX, DY, XK, YK, XKM1, YKM1 
CHARACTER*1 MATR(53, 105) 


x 
~ 
mot 
»x< 
o~ 
~A 
Ww 


XKM1 = X(K-1) 
YKM1 = Y(K-1) 


CALL COOR(XK, YK, KX, KY) 
CALL COOR(XKM1, YKM1, KXM1, KYM1) 


}ABS(KX - KXM1) - 1 

1ABS(KY - KYM1) - 1 

PMGi. Lk ) i = 

XK - XKM1)/(L+1) 

YK - YKM1)/(Lt1) 

= 1,L 

= XKM1 + J*DX 

= YKM1 + J*DY 
CALL COOR(U,V,1U,1V) 

1 MATR(IV, 1U) = ‘t’ 


OS 5. 
SLX 7 
Wot 
= Wusa 
rN TN 


Se 
Cc 


RETURN 
END 


SUBROUTINE GRAFOB (X, Y, NOBS, MATR) 
C KKKKKKKKKHKKK KKKKK KKKK KK KKKKKK KKK KKK KKKKKHK KK KKK KKK KK KK KKK HK KK KK KK KKK KKK 


REAL*8 X(1000), Y( 1000) 
CHARACTER*¥1 MATR(53, 105), CH 


MATR(27, 53) = ‘*' 


DO 2 K = 1,NOBS 
XK = X(K) 
ie ck) 


CALL COOR(XK, YK, IX, IY) 
CALL CHAR(K, CH) 
2 MATR(IY,1X) = CH 


RETURN 
END 
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SUBROUTINE CHAR(K, CH) 
CRE I ER RE RE EE EE A A OE OEE OEE 2 2 IE OE I 2 OE I OE I I A IE OC IE CE I CI 2 2 IE 2 a oe 2 oc 


CHARACTER*1 CH 


IF(K.EQ.1) THEN 
hie ali: 
RETURN 

END IF 


IF(K.EQ.2) THEN 
CH = ‘2° 
RETURN 

END IF 


IF(K.EQ.3) THEN 
CH = 3° 
RETURN 

ENDIF 


IF(K.EQ.4) THEN 
CH = ‘4’ 
RETURN 

ENDIF 


IF(K.EQ.5) 
CH = ' 
RETURN 

ENDIF 


THEN 
5' 


IF(K.EQ.6) THEN 
CH = '6' 
RETURN 

END IF 


Ie CK EO. 7 ) 
CH = 
RETURN 

END IF 


THEN 
a ' 


IF(K.EQ.8) 
CH = ¢ 
RETURN 

ENDIF 


THEN 
8' 


IF(K.EQ.9) THEN 
CH = '9' 
RETURN 

ENDIF 


IF {KEG 10) HEN 
CH = '@' 
RETURN 
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ENDIF 


NECK. EQS 11) at HEN 
CH = ‘#' 
RETURN 

END IF 


IF(K.EQ.12) THEN 
CH = '$' 
RETURN 

END IF 


IF(K.EQ.13) THEN 
CH = ‘%' 
RETURN 

ENDIF 


IF(K.EQ.14) THEN 
Ch =e 
RETURN 

ENDIF 


IF(K.EQ.15) THEN 
CH = ‘¢' 
RETURN 

ENDIF 


IF(K.EQ.16) THEN 
Coane, 
RETURN 

END IF 


iF (CK EGE 1iFEN 
Coa 
RETURN 

END IF 


IF(K.EQ.18) THEN 
Chie 
RETURN 

ENDIF 


1F(K.EQ.19) THEN 
CH = ‘>’ 
RETURN 

ENDIF 


IF(K.EQ.20) THEN 
CHa s 
RETURN 

ENDIF 

END 
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APPENDIX B. EXAMPLE: GENERATING CLASS NAMES 


This appendix illustrates the procedures used in 
traversing a ring graph in breadth first order to produce the 
names of candidate homotopy classes. The traversal is done 
Eollowing the rules set forth in Chapter III. 

The problem considered is that of determining the names of 
all candidate classes of paths whose shortest representative 
does not cross itself between points a and z in the plane 
containing two fixed obstacles. Figure B.1 illustrates the 
reference frame for the region and the location of the points 


aand z. 





Figure B.1 Region with n = 2 
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Once the reference frame is established, the region is 
modeled with a ring graph as shown in Figure B.2. This graph 
is traversed beginning at node III to produce a list of all 
character strings representing class names and the vertex jan 


which a path in that class ends. 





dBi | 


Figure B.2 Ring Craphetor m-a7 


For illustratien purposes, as the algorithm 1s 
implemented, all k+l long offspring of a given k long string 
will be listed. If a k+l long string violates any ot@eie 
"“sieving' criteria presented in Chapter III, excluding 
consideration of the ending node, it will be marked witha (*) 
and the rule it violates will be stated. That string “wan 
then be withdrawn from further consideration in the algorithm. 
In addition to the criteria set forth in Chapter III, those 
strings which duplicate any previous string will also be 


labelled and discarded. 
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Upon termination of the algorithm, those strings not 
labelled with a (*) will be searched for strings which do not 
terminate at the desired node, node I. These will be labelled 
Veen an (n). Those strings not labelled with a (*) or (n) 
upon completion of this final check will denote the desired 
list of class names. 


The list generated by the algorithm follows. 


Verolat1 on End 
Code Node SEming Violation 
tad e 
it Qs 
ry a) 
i Oo 13) 
I O90} 
* I Oj Oy duplicates aa, 
I a) 3, 
Tel yf) 
lay A523) B, 
nh IV 913) > 
n Id OO, 3, 
* IV OQ) > cancellation 
n IV Aa) By 
i IV O)B4 a> 
n it 1 Ba) 
n tok a) 8,8, 
* Te OB) a) B) wraps 
n lek 13) A) Qo 
* I OB) Boas wraps 
n ri] 03) Boa, 
i 913) a» 24 
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From the above list, 


class names are: 


1 & W NN FF 


OB) 
age 
a )B, 
0783) a8, 
a) By a)B, 


ee 


Lit 
Tt 


Era 


IGA 


eel 
IV 
it 
IV 
ie 
Pt 
Ta 
ie 
IV 
IV 


Piet 
ei 


OB) & ray 
O70) Bay 
O70) By» 
OQ) BQ) 
OQ) Boas 
0) By O90) 
0) By a) B» 
0) By a) B; 
0B) Oa 
0) By Bay 
0) By Ba» 
O93) Oy Or) 
913) Boa) @ 
0913) O93 a» 
0913) 939 a) 
O73) a) By By 
Or) By) Oo 
Q) By @)B) a) 


(0) B5a)B) a, 


O39 &) 1B) Bo 

0913) & 9B) Ba, 
O98) &)B) By a, 
Q) By &)B) Boa) 
0) Ba) B, Ba, 
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duplicates nf3) XO, 
wraps 


duplicates Oy f3) Aja, 


wraps 
duplicates OQ) Bra) 


wraps 


duplicates OQ) Bray 
wraps 

duplicates OB) Ba, 
cancellation 
wraps 

wraps 


crosses 


wraps 
wraps 


crosses 


crosses 
wraps 
wraps 


crosses 


the remaining candidate homotopy 


Figures B.3 through B.7 illustrate paths which belong to 


each of these classes. 
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APPENDIX C. EXAMPLE: COUNTING NON-ISOMORPHIC GRAPHS 


The material of this appendix is presented to illustiage 
the ideas introduced in Chapter V. The example used counts 
the number of non-isomorphic ring graphs given n=7 using 
Burnside's Lemma and then demonstrates the procedures used in 
determining the number‘of ring graphs which may be generated 
given a particular partitwon ef nme 

Given the number of obstacles n= 7, a regular 14-gon, as 
shown in Figure C.1 is produced. Edges are numbered clockwise 


for convenience. 





Figure C.1l 14-gon resulting trom n-7 


The edges of the polygon are then colored black (bold) or 
white (fine) subject to the restriction that opposite edges 
are colored differently. The question to be answered is this: 


How many non-equivalent two-colorings exist? 
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The 4n = 28 total permissible permutations consist of 
2n = 14 rotations of the polygon, n = 7 reflections about 
diagonals drawn through opposite vertices, and n = 7 
reflections about bisectors drawn through opposite edges. The 
cycle notation for each of the 28 permutations is shown below 
along with the number of colorings left fixed by that 
permutation. These numbers are obtained by applying the 
formulas introduced in Chapter V. In cases where no coloring 
is left fixed by a permutation with more than one cycle, 


cycles containing opposite edges will be underlined. 


Rotations 

(1)(2)(3)(4)(5) (6) (7) (8) (9)(10)(11)(12)(13)(14) = F¥(m) = 2’ 
(2°34 5 697 8 9°10 1) 12°13 14) = Ft) = 0 
(1357911 13)(2 4681012 14) » ¥(x) = 2 
(147101325 8111436912) = ¥(xn) = 0 

(195-9 13e3°7 11) (2 Galo 14 4.8812) -2F Cry), = 2 
(161127123813 49145 10) = ¥(x) = 0 

Ge eeSell 3 SiO.8e14 6 1264810) -4Fo)e=02 

(1 8)(2 9)(3 10)(4 11)(5 12)(6 13)(7 14) = ¥(n,) = 0 

(109 3 TieS2ls 7)(2 40 412 6 148) = ¥ (7) = 2 


(lO 5 W4n9e4 13°853 129792 11,6) = Peg) = 0 


SS Meee Gee creat ee cere 


mM : (11173139 5)(2 12 8 414 10 6) + ¥(m,) = 2 
yee: (1 12°9 6 SelMpeesa2e13Ml0n7N4) ake) = 0 
Mm, : (113 1197 5 3)(2 1412 10 8 6 4) = ¥(m,) = 2 


mM, : (1141312111098765 43 2) = ¥(n,) = 0 
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Reflections About Diagonals 


m, : (1 14)(2 13)(3 12)(4 11)(5 10)(6 9)(7 8) = F(m,) = 0 
TMs : (1 2)(3 14)(4 13)(5 12)(6 11)(7 10)(8 9) = F(m,) = 0 
Mm, : (2 3)(1 4)(5 14)(6 13)(7 12)(8 11)(9 10) = ¥(n,) = 0 
m, + (374)(2 5)(2"6)C7_14)(8 13)(9 12)C10 11) = Fea = 
Mm, 2 (4 5)(3 6)(2 7) BS 14)G0 13) 129 Ya = 0 
Mm, : (5 6)(4 7)(3 8)(2 9)(1 10)(11 14)(12 13) = ¥(n,) = 0 
Ty : (6 7)(5 8)(4 9)(3 10)(2 11)(1 12)(13 14) = F(n,) = 0 


Reflections About Bisectors 

wet (1)(8)(2mi4) G 13) (4 12) 11) (6 lenie 9) = ia 
Mh, : (2)(9)(1 3)(4 14)(5 13)(6 12)(7 11)(8 10) = ¥F(m,) = 2 
%, : (3)(10)(2 CQ 5) 14) 13) (62) 11) = Fie = 
Mm, : (4)(11)(3 5)(2 6)(1 7)(8 14)(9 13)(10 12) = F(m,) = 2 
Mm, : (5)(12)(4 6)(3 7)(2 8){1 9)(10 14)(11 13) = ¥(x,) = 2 
Nm, : (6)(13)(5 7)(4 8)(3 9)(2 10)(1 11)(12 14) = ¥(n,) = 2° 


M%, : (7)(14)(6 8)(5 9)(4 10)(3 11)(2 12)(1 13) = F¥(n,) = 2 


Applying Burnside's Lemma to the above permutations gives 
N = (1/28 ) ( 2' +2 4a ora ee eae 
or, N= 9. 
Thus, the number of non-equivalent edge two-colorings of 
the polygon, and hence the number of non-isomorphic ring 


graphs for n = 7 fixed obstacles is obtained. 


Next, we consider the partitions of the integer 7. They 
are given below with an arrow (~) indicating those with an odd 


number of summands. 


~ 7 
6 1 
2 
- 511 
4 3 
-~ 421 
4111 
-~ 331 
= 3 22 : 
Beek, 1 
eeeo i 2) 1s) 
222? 1 
ee 2. 1 
Pa le, te) 
el Le ll hve 
Consider the partition 2,2,1,1,1 with k = 5 summands. 


Applying the information from Chapter V, a pentagon (since 
k = 5) is constructed as shown in Figure C.2. Vertices are 


numbered clockwise for convenience. 


4 5 


Figure C.2 Pentagon Corresponding to Partition where k = 5 


de 


The 2k = 10 permutations of the pentagon consist of five 
rotations and five reflections about axes drawn through a 
vertex and bisecting its opposite edge. Each permutation is 
listed below, along with its cycle notation and cycle 


structure representation. 


Rotations 


Mm + (1)(2)(3)(4)(5) = x! 


Ty tel een ede axe 

iy “Cle Sw on2 4 mera 

mc : (12 eSes) ee 

Toe dS a 

Reflections 

Mm : (1)(2 5)(3 4) ~ x) x, 
(153) (2045S) =a; 

Ty 2 (1S) CsA sees. 

m% : (1 2)(3 5)(4) = x) x; 

My, 2) (ANIC ZS (Elec 


These permutations result in the cycle index 


Po = A(x? + 4x5 + 5xq%Q) . 


In applying Polya's Theorem to determine the number of 


non-isomorphic ring graphs generated by this partition, the 


2 


number of colors used is two, since the partition consists 
only of the summands 2 and 1. Letting these colors be b and 
w for black and white respectively, 

Po((btw), (bi+w’),..., (b’+w)) 


yields 


= [ (b+w)® + 4(b5+w>) + 5 (btw) (b2+w*) 2] 


From this index, the coefficient of bw? 1s desired since the 
summand 2 appears twice in the partition and the summand 1 
appears three times. 


Expanding the index gives 


= [10b° + 10b*w + 20b°w? + 20b7w° + 10bw‘ + 10w*) 


Hence, the coefficient of bw is 2, which is the number of 
non-equivalent two-colorings which may be produced from the 
paGertron 272,1,1,l1. The two corresponding non-isomorphic 
ring graphs are shown in Figure C.3. This same procedure is 
applied to the remaining candidate partitions to find the 
number of non-isomorphic ring graphs which may be produced 


from each one. 


13 


Figure C.3 Non-isomorphic Ring Graphs Produced Using 
Summands 2,2,1,1,1 


114 


LIST OF REFERENCES 


Willard, S., General Topology, Addison-Wesley Publishing 
Company, 1970. 


Thornton, J.R., Algebraic Names for Homotopy Classes of 
Paths in a Plane With Obstacles: A Foundation for the 
Shortest Path Problem, NPS Technical Report, Unpublished. 


Cuerington, A.M., The Shortest Path Problem in the Plane 
With Obstacles: Bounds on Path Lengths and Shortest Paths 
Within Homotopy Classes, Master's Thesis, Naval 
Postgraduate School, Monterey, California, June 1991. 


Carrol, J., and Long, D., Theory of Finite Automata, 
Prentice-Hall, Inc., 1989. 


Dossey, J.A., and others, Discrete Mathematics, Scott, 
Foresman and Company, 1987. 


Buckley, F., and Harary, F., Distance in Graphs, Addison- 
Wesley Publishing Company, 1990. 


Tucker, A., Applied Combinatorics, 2d ed., John Wiley & 
Sons, Inc., 1984. 


Hillman, A.P., and Alexandersen, G.L., A First 
Undergraduate Course in Abstract Algebra, 4th ed., 
Wadsworth Publishing Company, 1988. 


Fredricksen, H.M., and Kessler, I.J., “An Algorithm For 


Generating Necklaces of Beads in Two Colors," Discrete 
Mathematics, v. 61, pp. 181-188, 1986. 


mS 


BIBLIOGRAPHY 


Patterson, E.M., Topology, 2d ed., Interscience Publishers 
InGc eee 959" 
Read, R.C., "Polya's Theorem and its Progeny,” Mathematics 


Magazine, v. 60, no. 5, pp. 275-282, December 1987. 


Roberts, F.S., Applied Combinatorics, Prentice-Hall, Inc., 
1984. 
Simmons, G.F., Introduction to Topology and Modern Analysis, 


McGraw-Hill Book Company, Inc., 1963. 


net -6 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22304-6145 


biorary, Coge oz 
Naval Postgraduate School 
Monterey, California 93943-5002 


Commandant of the Marine Corps 
Code TE 06 

Headquarters, U.S. Marine Corps 
Washington, D.C. 20380-0001 


Professor John Thornton, Code MA Th 
Department of Mathematics 

Naval Postgraduate School 

Monterey, California 93943-5000 


Chairman, Code MA 

Department of Mathematics 

Naval Postgraduate School 
Monterey, California 93943-5000 


Professor Kim Query, Code MA Hk 
Department of Mathematics 
Naval Postgraduate School 
Monterey, California 93943-5000 


Captain Kevin D. Jenkins, USMC 
712 Lisburn Road 
Camp Hill, Pennsylvania 17011 


Captain André M. Cuerington, USA 
P.O. Box 114 
Highland Falls, New York 10928 


Commandant of the Marine Corps 
Code MA 

Headquarters, U.S. Marine Corps 
Washington, D.C. 20380-0001 
Attn: Captain D.W. Hintze 


7 














ha Ui 





Thesis 

J4565 Jenkins 

Ge. The shortest path pro- 
blem in the plane with 
obstacles. 












Lat very ry i ee Cees 
, ’ re, per ai iS ety Se renee Ad A Beat, t ] 
he 


labia ee Cou ate SCC AEE Pe LT Ee fe 













an yt a eaetag 6 
4 fe J Fs 
rene oes aes ao Ra 


ea re eee) ees . Ppa, + - 
af aoe Se aera a ae te eee 
pie 





¥ ca Sane He DU a 
Digighservlot ; es nee es a Re ND a : oN) ny ; 1tgie ay al iH 
a QUdY ci te Nard Sat: A NYY cyan more pure 4 1 fl il TANT 
eer 7 2 : ev 9 nee Bat ae eT Fm eae fe Th | | 
baccarat eras meh id. Mids coikar ey a ameay Raat He TUE il MI 
pee are a f cares aera TIBI S EER hy abt shi rset oh aa : 

a trdh Semen ora Gateceressh co ne o ey rc ee , es ah tow: mK tt WW | All ith WI} \ Hulk 3 | 
ARE Jootar ron ne OR RN are Aes Prope cotree rg INE yee . 
Peet treet ee Sen Hager LK APE DE Ras Agee) coeur a re Oe Seas a 3 2/768 HAY 6358 
a Be kgs tet ary — 


ota 





























































































































aie 5 a) Ties aA “. ah Dy Sn . 
ne ae eran ov brea ee righ pea peat rae ry OE £ gamte! f r 
if fer Rect pe TONEY ia a See ot at Tae nn: ra 
rae Ronee ae oe SF ee Ras fie ere $a ves “i * ys De lta bees 
pines cra a ’- Rae ear Ary Se Ro woh eee ote nie eet oy. y AYP ; ; 
Pte ao one eH ieieeen retary AE aC SECO ES ed ode OES Se Ake - 
eer sees a ata pias RT Ce NaS “ SUA ae nr ee bis. ’ , : . 
rr ar fees yori alt oa & prety Te Pee, eS ; 3? 
ord ye po ots Bese Ci ic are naa sore ahs Beer Leer Ot Ao wy "ew F D m4 rs F U D e 
SR er Co wenn Patt RTH We} ek 3 wy Wiss ’ ‘4 4 ay 
4s rr Rts ree} +d wen ea : eae Rees Rar ee oP Pe ae mi Ge A é x Fr : 
ne aoa ie 9 haat areas) ger are Bigeye doh iS 5) a get ee wo i q J vu Fj 
} Prarn ree ee ia Wa Ore Se PS eG oy nes Reon Ce wn Sad Leet — A A a = F 
? mieten ito eee ae SF Sr AREA o} nob ys be sta EE Lelia ar: +t, ‘ ae & fh) ae - 
Petian sii eoat sites Beri Les eae ae OY rea masts Fi veers ie : age 
Serer) - ga er ee te tae ype . 8 5 : : 
ae ota ee RCN DATA ate Lee . # i eens, Ms eg lae y f OE 
ty PVs te a shades pee $ tf. “e : . hi 
ah ee En ks ae pone Ces Pa Te a ne ta : C p *. F ’ ar ‘ F 
ape Se pate A Prieta! Get, ee Bh ts bs ' - a 5 
ees ri cipey, tae veri ro ee, Ne " eet y . a 
Ne i Be ve. sat cea = : 
hh oblr runes) ve avs Creare zy bed ne rhe is ah > is rs bY 4 : 
SY Decca rae : CO sean tren oye a ) F 2 : 
andehe tise. a chet Shea} ee Ress iure es Te he eS eu ete 5 . i oo 5 + P 
se tee a shale Aeeneee ae = Bork: Ala Sea ee et ne brad be ar gly: % a a 
D ey v rare. E une. oS) rae fl bell" Poses ate Bn ae BR eet * def le ae . soo? 4, rr ; \ : 
1 ae i z rr sera ee PE at ae 5 Pare ae 
1} 5 pein et seanere Ue tr ae Gragh od eM a Ae Pe cS ne aS) : 
pays pA re pare AS ee Pipa gh te Ay by p ie cae ra eT Cer tA ona rae Ae es ea ae ae zd ee ya i : " 
Sens Poh . are eae rte Ary Bott ed op Cy Co oe Fo et as an ae Wd . bd 7 
rr he as aire m ae a oo pat Ma 28 Seb tS a Pe ce ast Re MY “the : q Pa : : 
Pee a je deter oe Hosts 747 cree Teas Ay Te eT » eo ; 5 ; 
we Sr cae ny by oo Ages } ste DY ney ha fad abe ddientt) hie eS ' a x - ren D 7 
et CR eae ae ae Pt = Ce tye care] Chiltan Tite ee) 7 ae « 5 » y a 
pt TST Sagal LEER a aay tHe rst. : A < fi : : 
Sorter at IAT ah Se we ear oe Sule aro De Re ed ora 7 ae ; A da 
aro or feet ee TA os Pe pag Car haa Te PNY | A ihe ‘ an F : ; 
Cee aie Se oe eee hein gee el esis MRE a ns ae 
i Eel chad ram 8 4 te Pe ‘4 oe 4 rot Z 2 J ‘ 
, whe 2B Sy Saree ar nes ts ne Perey os i AD He ny OS t D : 4 * 
La J para oe ARES Fone ah rs Sia RA AES aca Se 4 oa a4 aah 4 He HS 5 n P yi P 
i pee te Ws aa i” era aS ), in ho ee Tie ' an bs v. 
et cee Tare! price ates D Bee ers nee oe * P eit , 4 & ; t D bs 4 
oes sans Fae ma Oi os te BAT aD fo en as) vas Te rr] } sims teks fe 
rere es soe cd ey ts Ras eS 4 ; ae PRA oh ¢ t F : RS ete 
oa 3 ra X ye . aT =e 1s : 
-< Saar eniae re ad ae bs 1 %&S % Tat Si p Sep” a me “ts r v By ] sg Sey m a . : 
Son eae Parte ss: aoe A ein eg 7 . J 4 + , ; 
<A a tle Baas <a ies nee ong Ada pede as in 5B 2° 
epee Baia ad Re Row rea 7 re € A 
5 wo shes re ; ie (ote nlf , age -— 7 oa D o 
. a vi “ p i 4 : y , F « 
Ps a Py . A fn 











Lh ree 
ne 7? me $5 
at 


ati 
aed . Pig tae 4 th ons te A 
Ses Petree eh ap dope n A tn eee 
ee ts 9 ic SMe rbey Sh ae DS dade ee od Etat ent and 
este Se oad Or a a he] Y 


4 nae rs ehrit 4 CLT, 5 
ig ae Bot ‘ Faet ane a 3 coe a ay ri 


























rt 
eS hein es ay 
ES Te 







Ha irk i ee Ee 


1 eet ors 


Ie pas, yea eart a ars 
oe re Pato tg * rr tap Tee oe eee hes + vf eae 7 
ot PTs % $4507 .P %, 
Peas a fe ao sneer nag vs TAKS sinc aS As 
Seca esata ee Ree eee Pa Pa CRAP LE a | ) 
2 A di - ‘ ‘ad wef S oo a - A - f J “ : 
re ae ieee are k, mh pee eet ore inte ak A Hits ig poe 3 ere eee ‘ Pes lm ae ; : Pare it 7 A 
io Bh figs Perr E + Butane. y ; 4 
sean sitet oF tears oy 
Porn a pe ee vie a ifr ee Teg oP) 
epee iie a cen te ae 















See 




















es ae eed 21 TS 











¥ a 5 
cep 0/ i 
i Rees br fy E 
“s Tek pia ey 





mE a) 
aed 


ay an Ps eh 
ie Bis ie ie cs Rites 












4 a : i es a Y 
Pe ole aT th ok hal de 
Pitt ems 3 ai ail 

a 4 1 2s ‘he, 

Wier hate 84 Po 
aay Se 408 cd HE Ag Bar a me a 
pete 7 We es oie ar arate ae Sy. gee ey Bey ft? 


iA = 2 SN : int be 

Riegel un >A ee 
bef Phact fx, Baad ed ' i 7 

md tbe erry saphena y, ae ie Y seat ed Rn eae ee a oa 4 
gay in Ata! Sty Crees eg Ee a eae . ako yy a a 


ie) > EAE RTs St) 
DRIER ret Meee o Shia aig pees 
a 
4 
paar S iat 7S paren uA a: ee : Baie oie mona ETE ee aI 
eS rhe te aes 2 ai ras reads sla ane fens CUP PLY oe Pt) oe ity Pik ataite Wad ToT Be 
Ee ot ert cee eh Sat ahh ae he 4S 9 uy b ee fe) Lok 3 a AS “oe <4 ae 5 
Bieri aie neere STS Dei PaNe et oe ae ve st HR Se ae Ry 
a pegs ard ey ae {te Rehr ae ag apc ve be Sa ee et he 
os Rage Tras Ore ae AL aa cases 4 7) ar: 7 aA Les at a + 
ted ¥ beh od Ne fied Pu MN Lay c ee ree Perf an 3 Net 4 Z RL se te fog SET it by by 
apt tg ah la nin Pana Renn A ea ee ences en, tie 
| cisco aul tit ae wrt 7 parent oe ; Pater take mre ‘gore irks +38 SOAS To 
Ag 3: o ee Sek: are 3 at ‘2 » ‘a: ab })2 i 
pee ae eee: Be tid d Pac +i ¥. 



















































i } 
t Conia Si doa’ fi: Jos a%pi Ae 



















Ey} va LP ae ee Phy 
elias rs peas ata ‘ a4 z ie Des eat 
aaa RES ey a e +] LAL be ty Se 
aS ei rity ae x va Sept) aT ee at 
AG SiS Foy os sayy Le a ea hs A ie dn 
hoe are Ret Sg ae pe 9 . NN + ee rea etna. 
oe om ee Any Ree , et fetiege fps rah cent % A es aS Liter ne " 
Fae aeoa ee ron Se Ta ie ee ae Be ein as att UO Shae i a A Ere ae Rest b has As hg 
he Org a , Sek tb! a kde aie 14 ded dG as* Cee: 
ve b ‘ipeinie R E Ako es ak a a ee tos Beet . ee Ta : 
Se Ee PR Rete) i fy PAS WET CNET Ne, fee Ver puss! ats ayo 
a ae te Banke cleraey DU chades sa tas Great Gite ae me Bhat a at sett ths oa 
Lae oe Lae tas ate > aD hed cee Pe Sc 3 Pars LC a 4, 
SeSter Gale ie eH rs fea Le res = hn > 
tee nt fe ape ltD TE re ee Ses 
. ov ar 
Yap Poe ee aulp od aatd se bate 
Rene ta aoe 
rene Sia ate 
: Se ere rg ean seer 
Red ky BUN Lin) 
¥ 
ng paula eer rs hoa 
Pears Mears 6 ata se 
behe ta aetna ak reaver iT “ath Fee ee 
rae ed a eae ceed aa 
.¥ Qt 
*: lads oly f-7 e* 


a ik Sy vi ety ve 
rin taal Te lnc tes ate aceasta: ie iG 
nae toreeiet gs AOR S 8 atc eed aise a 


Ne ate ay 
tre 7 find f : 
pone thd dhe opto eamA Ly eta ded cele gd peepr te ait pte - 
nee gree RAE tir ie he ee F at hte ened Sta (es a er) ata a 
het ae a hd Wen acs sat od a + Pars cal eae ig ag hy tents bh wae 
eg tae ai aa Sa Sbeste let < toes Rstrnt aire Laan Ne Te De Soe Tt cho ne 
ie eS eS eS At habe see} > 
= ED Pedro Saya ba -¢44 ee ak sot eee Fe Labs bk Tesla ee Ph 
Ee Negegey Soret eats & aa ae Ret ts boety iva ote dt ete ata a 
“HP red a oa as] es ata PPTs. 2 oC Ps 


agttecad ite, rr EN pi in iT eS pw ha a Lai ich lhe AL debt 
Fe gree Ue ge Ran pee RN 
= foro inet oes a Uae the 9 ae 4 bin bay 4 Raa thts 


ene th ae rs aed Seer ngh: peer ffs Behe Capt a ve te) ibe Pe 
ef 5 i oA AT & Laie Dieta teed deed 343 hd 
ea Nia getdate de bank aed Wd ia isetead Ms Emcee atl Co: Red eee hi bth th We Aled 
Th ah plat “bee Cok, littl ilk aed ahora teat [epee eee ae td Pane tale od Oe wipie ey 
Wins 1 ea) nak Tal if writ hh ten A parse} Y 
a ; ET ia ad Dat elon fe LS et Aas Cae ef soc Ne 
Tt SN vcs Soy tases fete Pee Le it del ro) re 
Pa aa od med es 5 Pt 
eRe “ca ath se Sie pa ote td Lrek une mess tie A ae a 
¥ Fama : J . rr) eas J » a a 
Feet era aaa ie pe ee BAR reat eh nee 
Seti aeetes Fa: > » PR “ga ta aha) 
; ice “* porate Sera eetfat bikin ve ay eta Prat EOS 7 haa 
Seater bake fd ies iH = 4 = 6 
Bota i Soret Piper Scted ro land dria bard Lia Sey ap ee ag 
eee Sher fn Pi bass Ad Net led Fad 0 cod jus allel top, net id, Cede ed ) 
nes ae ae e, : 
hh ae vf ah pts ie 1 fy Feel ied AF Gd A he 
es aeeese sebetrens Cena ae anit le 





Ro ; ise * 
Aiecks thts ge S 


















Bos teh “4 
eer nb ay Pr : 
ah a Rete a 
Pee Drs se ye te ar 
Rea aco o¥e)) = 
Soret Ra SOO eee te el ae Ok 
wey betes a ee A ff. 
“F 32 q'ee9} 
Ltad 
























od 
7 

bi 

4 

re 

rd 

ry 
xs 

a 


.* 
























nb 5 


rh 


af 
s¢ 4 Ca ah he Lael Od de pi 
es ia 


oe 








pots wren 9 2) yeas aa 
Pe tek ol aE, ax Oe oe. ae 
pie ae eo ‘s ager Uy ew Alco ry, ar i 
a are Sc onc a sha tres ae aL sp apse es GM Ta he ls 
EY chabalel Ty: POs ee ch rn fo PELSOR COMO NT OLure Beghm eee | aE, ee Je 
2 pes tats eee J stl ie peo aA aot ast papas Ea el epee 6 rat te tt aed art eS ~ i 
alll i Seda Sate d fale’ ¥ bh te 5) wrayer Le ed Me elated el eae ew - 
ae cennpay eee Lrg a Ree SS he Et de die an tone red 9 eR et oe ay ARES ® wae 
whe Ke erate pam PP tan Chi eon at PIE CAO, TP ages 7 i 
iad. Conny hae loth ld rien: a ad eet ae ee Ak ae a th 
a feet ad) ane or Fda us re Pre York ik on ae el be hae Ca ck hatin Bil 
Oa Lidantal fun aby atl E | oh naiiie Sie, eat Phe, ra ta a de St) ea 
bel epee pr ian cds yo Fahy be th sania pest, |Add ch de fag "ee wees \. 78 wre WTR: 
Phpecrestarae he nen, Pro tS, died eee Bas ak dd hich o hashed Chr oe tere are wh, my m7 
pataforeegincuabag rit phemeaeneeh Meee De etree tay oti ya ek aes be Do ae 
Fy et tae 6 Ph ered oe: Ne en ome oe tn Fete he eh ee Skee Me es saute ts rw - 
ne nannies foram. fe pe aie te +" Bh meade for, nA 4 i raps ny , - 













3 he sfeens ai aad ed a pa Pty oy ia ‘ 
ree pee de hla ela hp plies tc Papi Gag eis pe fee conte eM eee of SE ye i Ra 
LEO tpt ae hn tlhe? sa oes earentanal Foca t x etate ite be aie aka c aa 6 eR Lie ae Pet ery 
hd De tled Lt hid 3 Paes lag le Pad ea Cal ae Leh ous, v, Tatil FL’ hk 7 asckedilsdieitin! tila ed, ibd Des dal 

babe 2d Seta eke ks Soe iat SL. hi Aqua hats) Aah Aes LS, 
eatgetem Lipa Lied merpyrenes aR bie tf re <r oa Ala dS elle a ad Pay tae, Ti) oh af, 
ae tame Ld far hang Co lag er at eae ok Packet tae sg A eee 
‘ ese id ng oe yc bor ery x Petras a os merece 9 er ede ed OF cs ub hag ad 
Setcncmbes eet eres Hcl Steer G8 View weal ae bere Ci cai ie ees fod rs Bear WR COOL Pave elites 
SY ; oils 6 nets 5 ck LT os giveth eye 
oo : c 


Rs hla ed 
rane a tk or Fig 
ae Seen. eine eae ee fa ea eee St nals a Eieon heen $s pena 
1 lhe hig da tea ay Ls A . bg beta ] 
Ree TED el eaten eee tit cenane See eee mT hel 
oe Sn etal Da: Sto Wa he tes 1 nee sh es Px ho 
in do tolcs Auth ‘gp, 8 SORA fon cunt Fett tp yA reali BCs th dtd 
> Bt ons abated eget ane ha et) tl Tondo tah Pbldh’. dk doe Pe hae i ud a 
e aici ean cS MR ech autem PT hina Cele By Dick ger otek Fae 
ala beast epee? Pei Ns EEE Shas Hid ee piesk teak hh Teal ene 
3 PEt snetealens Sah ee eae pecee te See pre Pe, LM A diet Madbei tet meant f w 418 ks pea a Td we. ,¥ « as Veta 
AF nag hry ahaha ithe Sc gape Se aE AL oe belay erate. Ad Ct par 
peettieea ‘el is aa eae ele ee i edn 8 ig a LES ea eens i 
Adis " are Laat pet i Ea as Bh “aeptes oa Neg. eae Abana € vee Ahh on. 
eeep wy a are (a Ad y MS hE) hae hk PA tory B, 7 i 
nape pa ibe Se baled ta ib dbs y rT) ve bd Xk ta Dent dada etl AT ae ry reemay 78 
faa tata eng e ‘dnepuadep Regi falda} Da ryee yrs” a eae Di he / dh Ms hook hthed ta Pe LE eed by A} OS Fi teUS @ 
Gea? ademas Ai Met Ne hie nk Pad a lat a ee Oe ey eae * 7am hc 
ie Aap ede a EN ees EEO, RR TA 
on " far ys if o Jit a0 ewe IT OF Rats an Id " ¥ Le 
La tee a Rana SLOG epekaal mean wre Tene ace ale Ce ee TO ; 
6 ae PEaat Ra ane wis Betta ah eaey ve i i Ks os Th Labs ake SX Ve eat Re ee Le Bg ae 9 ; ’ ‘ 
et oie td shal oP acid CINE a ee to a aad no Pea a ATA eld vy VK ¥ LAT) aha) me 4 
; San bea ie PNG a habe neat mre eg s4 : Foe Tsk Ries A ‘=r |? *'. 
es bait Hy 5 Ie Lada ats eat eu “ wey * Al Seek ee ee D LP hy bs ‘ ® 
Lanne) ma A ae Cea ath eae ube ang eee 1 f me 2 4 


















a , es ” . 




















panics thee ey dhe Saari Le tknd 





Dee Saad 
ea Sree fa 


rete Na tae calla 4 
pita ea Eee fe 






Pe ee ke ee Bt ons he ve Pre iets as ae Qa ad owt at | 4 
Biron ee eee pe A tats ei ap trae th PP ath ane, Odes wr P 
a Waa were rat tied it Sag at Meek Comme ar ‘bee ge F 
” ALS A SE OTS 
Ce ees <o% Enchehe. a a 4 
cele aie heh een pee, 
3 ae | 
re SO hoes Spee hs Sate FAL By 
pe Ntyren r=: md PT a z id 
5s hea Oe ale ab \ ) a ae 
AS Tha ey BEA DAA ati Tia AR SL SE 





yo 


